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ABSTRACT 


Uniform and isotropic mathematical models of the expand- 
Hie ulliverse usually predict an initial singularity of infi- 
nite mass density and space curvature. To study possible 
mechanisms which would avoid the occurance of these singular-— 
HG LeES , non-uniform eeomolouitcais medeils based on Himereim is 
Preid equations are investigated in which random perturba- 
Mions Of long wave lengths are superimposed on the Robinson- 
feebter i mnetric ror the unperturbed models, Techniques of Iluad 
PaouLence  UNeOry. Used TO.describe random fields by a hier— 
meemy Ol Central moments of the random perturbations; are 
@oplied to describe the dynamics of these moments. For the 
eecoe Of small perturbations the hierarchy is truncated and 
Bolucvions are found. The solutions are either growing or de- 
wayines perturbations leading to R™ extra terms in the usual 
eosmological equations for the curvature radius R. The re- 
sult agrees with the small perturbation Fourier series expan- 
* $10n analysis which exists in the literature. Based on the 
upper limit of the anisotropy of the 3° K background radia- 
tion, the growing perturbation model predicevus: 4a max imum 
expansion even for k=0, Euclidean Spaces. The decaying per- 
turbation solutions give extra terms of the form WR with 
m>4 in the cosmological equations and indicate that the 
mechanism of long wave random perturbations may prevent the 


original singularity and make oscillatory models possible. 
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oo Ye ODUC ren 


Perhaps one of the most exciting discoveries of the 
twentieth Century thus far is Hubble's now famous discovery 
(1936) that the relative shift in wavelength an iene 
mecctralvemissions Of each star iS proportional to its 
astronomical distance L [1]. Modern cosmology generally 
mecognizes this phenomenon to be the consequence of an 
"expanding universe" whose dynamics are determined by the 
Einstein field equations [2, 21]. The concept of an expand- 
ing universe presents various problems to those who are 
mavolved with the systematic study of our universe as an 
entity. Because the term "expanding universe" implies that 
feel adistances between all points are on the average increas— 
meee tt 1S Natural to ask If our universe was ever in 4 
fendition in which all lengths were zero. it is also natural 
momask if this so-called expansion is homogeneous and iso- 
treopic at all points in the universe. 

ieuene answer to tne farst quesvion 18 yes, Unen it is 
eeearent that sometime in the distant past the mass density 
of the universe was infinite. But Wheeler [20] has pointed 
Sueethat such a state is rather dubious. He calls attention 
memthe fact that it is possible to construct a "length" L* 
bya combination of the gravitational constant G, the 


light velocity c, and the Planck constant fh such that 





= bai = 1.6 x 10” en. Gel) 
e3 
memSUPBeStTS that this length might be a characteristic 

length that is associated with a quantum theory of gravita- 
tron. He points out that because of the smallness of this 
dength, quantum gravitational effects, though negligible 

mow, cCOuld in principle be important during the very earli- 
est stages of the expansion of our universe [5]. Hence, 
iets believed that the undesirable singularity associated 
reecir the iniinite Tenet ey MDigoeeOCraVvOl1ded. = One sea com- 
Memsays, why such quantum effects might be significant 

meine Lhe highly contracted phase of our universe is that 
Moe Cstimated cross section for pairs of gravitons to produce 
meirs of electrons and pairs of photons is comparable when 
Eme Wavelength of the compacted gravitational radiation, 
more is assumed to be present, is of the order of L*. In 
Sager FO GeSscribe such @ situation, it is clear that a 
quantum theory must be used to treat this obviously non- 
classical phenomenon [20]. 

Various workers have begun the task of incorporating 
m@eancum effects into gravitational theory. some investiga- 
tors (Misner [16] and Parker [17] for example) do this by 
werdous canonical procedures, but others have taken a differ- 
ent approach. Woehler [34] suggests that it might be fruit- 
mr FO postpone a correct quantum theory of gravitation and 
try a semiclassical approach. A way he suggests to do this 


immco presume that the field variables Ey (called the 





metric tensor) of Einstein's field equations randomly fluc- 
tuate with position and are characterized by a "mixing 
length" L, which for quantum fluctuations would take the 
value L*. By considering the average of the field equa- 
tions, he suggests that it might be possible to avoid singu- 
iteesOolucitons if the extra source term, whieh appears as a 
consequence of the averaging and the nonlinearity of the 
field equations, is dominant during the early phase of high 
contraction of our universe. Ginzburg and other co-workers 
[7] use a similar approach. They attempt to estimate the 
memecvional form of this extra source term by expanding 
Chis source term in a series expansion of the average curva- 
murme tensor. It is interesting to note that both the tech- 
Pemees Of Woehler and Ginzburg are reminiscent of an earlier 
technique used by Welton [33] to estimate the Lamb Shift 
Pyea semiclassical model that also involves fluctuations, 
Smee or the position variable of the electron. 

However, the point in the above discussion that is 
meemene to the introduction of this thesis is that there 
eme  voday serious considerations of using a randomly fluctu- 
ame metric tensor et as a prelude to developing a realis- 
meemaguantum theory of gravitation. 

ibmere is also another aspect of cosmology that requires 
m@remuse Of a nonuniform fluctuating metric tensor eu 
mes Consideration is brought out by the second question 
that asks if the expansion is everywhere homogeneous and 


[ieeeropic. it is a matter of history that one of the first 


LO 





MyaecMmOvouUO Gescribe Re Universe [rom ao Cesmoleopicalk point 
of view assumes that all matter and all space are homogene- 
eis ald  lsSOUropie On a sSuiriciently large scale. This 
Mieorevical assertion 1s called the Cosmological Principle 
ea@dewas tormulated in 1935 [2, 21, 22]. 

Although the use of the.Cosmological Principle appears 
Power lyoin tne liveracure, 1ts craticism is still very 
much alive today. De Vaucouleurs [31], in a recent paper, 
Ilernts out that because current observations have as yet 
mea led to demonstrate the existence of that sufficiently 
Wemee scale beyond which space is homogeneous and isotropic, 
it seems unlikely that it would be easily found (if such a 
length scale exists at all) by a slightly deeper astronomical 
Save y . 

Because Of this and other arguments, many workers have 
Paeweled TO consider memuniform cosmological models. ~ One 
Wem tiis 18 done is To consider the secs Of “an. anisovrop— 
meruniverse which has two equivalent and one inequivalent 
Garecctions in space. The dynamics of such models can be 
wevermined exactly by solving the Spoon ee field equations 
[27]. Another method to construct nonuniform cosmological 
models is to consider small linear perturbations away from a 
seavially homogeneous and isotropic cosmological metric. 
Sachs and Wolfe [24] do this and are able to express these 
MBetrlLurbations as a function of time. Their solutions contain 
meme arbitrary functions of the three spatial coordinates 


moreh are determined by initial conditions of a given 


pal 





arbitrary perturbation field. These functions then are 
mBenerating functions which correspond to a given nonuniform 
meael. 

Meat pae amportant polne for the reader of this™@= paper 
BO note is that cosmological models with a variety of non- 
uniformities are being actively considered by modern cosmol~ 
Geisus. ihnese nonuniformities may be zero point quantum 
mechanical fluctuations or they may be fluctuations due 
to a large scale uneven mass Gus tri Duy Olnes En racy wade 
beee and others [7] point out that both types of fluctua- 
meno are described by the same formalism, a metric tensor 
that fluctuates randomly about its average value. Matzner 
[14] points out that this formalism is appropriate for a 
Universe whose average background is uniform but contains 
Zewevational radiation of an arbitrary size L. . Clearly, 
some techniques are needed which describe these random per- 
imebations of arbitrary magnitude and size L in some natural 


way. 


ve 





TiO LEV PN Ge a bOSOni eon Dae MODEL 


In 1934, Milne and McCrea chose to describe our universe 
in terms of a Newtonian model. This model treats the uni- 
memse 2S 2 Streaming, uniform fluid. . Surprisingly, they 
emowed that such a model is in many respects equivalent To 
the models of relativistic cosmology [2]. Their model has 
Mme  adiStinect advantage of being mathematically much less 
eemplex than models developed within the framework of gener- 
mempeclavivaty. Another advantage is that iluid flows have 
meen well studied and if would seem, therefore, GaacGe Some wor 
miomoecwy tC SeOr 111d mechanics might be transierable Lo the 
meld Of cosmology. 

Such a possibility hase liepari. DCE Chemo iva Caen 
memepreparing this dissertation. In particular, the study 
of turbulent fluid flows which deals with random perturba- 
mens Of fluid flow fields would seem, on the basis of the 
above discussion, to be relevant to the study of random 
bercurbations of the "flow field" associated with the 
Maayerse. | 

It is well known that fluid flows obey the Navier- 
Stokes equation. This equation relates the velocity field 
V to a pressure field P associated with the flow. It is 
also well known that if both the velocity and pressure 
fields are decomposed into an average field plus a fluctua- 


meron tiield such that 


Is 





>> 
<V> + Vv 


<+ 
i 


rd 
il 


<P> + p (2-1) 


where < > denotes "average," then the Navier-Stokes equation 
ean be placed into a form which is symbolically identical 
with the Navier-Stokes equation for V and P but with v 
Replaced by <V> and P replaced by <P>. Owing to the non- 
linearity of the Navier-Stokes equation, additional source 
werms Which involve the correlations kl appear (where 

i and j denote the i and ee components -of =the fluctuation 
v) [4]. These extra terms are called the Reynold Stresses 
and may be interpreted as additional stresses which act 

upon <V>. 

There are a variety of methods that fluid dynamicists 
maecmoo GecLtermine the functional form of the Reynold Stresses. 
An important way is to generate a hierarchy of moment equa- 
tions in v by multiplying the unaveraged Navier-Stokes 
equation with all orders of v and averaging. Because of 
Mmaiesnonlinearity, each new equation so generated contains 
feeoment Of one order higher than what would be present 
if the Navier-Stokes equation were linear. Although any 
mamiber Of moment equations could in principle be generated, 
there always exists one more unknown than there are equations. 
hams problem is called the closure problem of turbulence 
@heory. Much research has gone into finding ways to replace 
this indeterminate, infinite set of equations by a plausible 
determinate, finite set so that useful information about 
furbulant flow patterns might be obtained [3]. However, in 


14 





spite Of the difficulties associated with the closure prob- 
tem, it should be pointed out that the turbulence approach, 
Moawehn 1S a sufficiently general formalism, can in principle 
M@eocribe any sort of random perturbation field of any magni- 
tude. (A short summary of turbulence theory of fluids is 
@ivyen in Appendix A.) 

The Einstein Field equations are likewise nonlinear 
Seer erential equations. Their field variables are the 
metric tensor ee and the energy-momentum tensor Tavs 
because certain aspects of the field equations of general 
Memacivitcy resemble those of fluid mechanics, it 1s nov 
Mieeasonable to presume that the Einstein field equations 
Seonla be expressed in terms of the average metric tensor 
Sor? the average energy-momentum tensor sae and some 
meomaual terms Chat appear in away analogous to the Reynold 
Pieesses Owing to nonlinearities. It might also be possible 
to generate a moment hierarchy which can be used for the 
Germermination of these extra terms if an appropriate closure 
@Meecmnique could be found. 

It is this speculation therefore, based on the similar- 
ities between fluid mechanics and general relativity, that 
feeares aS the motivation for constructing such a technique 
Women might be useful in solving perturbation problems of 


momoee relativity. Tne remainder of this dissertation is 


epee ply to this conjecture. 


ae) 





heel ria DeoOkI Ee? LOW OF aoe p iD 


Cavern Chav. Weeis pessible to penerate. a hierareny of 
Meaningful moment equations from the Einstein field equa- 
tions, it 1s necessary to truncate them in some manner. 

Then a solution for the "Reynold Stresses" can be determined 
faee Will be found to contain second order moments in the 
Womeurbed Pee and its derivatives. “Justvas in fluid 
Mmmgellence theory, the closure of the moment hierarchy is 
MemoriVvial because it requires special restrictive assump- 
mmenics COncerning the nature of the random perturbation 

field and does not just "fall out" of the theoretical fra 
framework. 

ine order to see how this truncation might be accomplished, 
momoider first the following three important truncation 
Emeeecaures which are used in fluid Turbulence theory to 
@mese its unlimited set of equations which contain moments 


> 
Gmeall orders of the velocity v [3]. 


1. Phenomenological approximations - These involve postu- 
lating (on dimensional grounds or by intuition) an 
extba relavion between various moments to bring the 
Mmeimower of equations up to the number of unknowns. 

2. Higher-order-moment-discard approximation - This tech- 
nique simply discards the higher order moments because 
of smallness or some other reason. This reduces the 


number of unknowns to the number of available equations. 


16 





3. Cumulant-discard approximation - This procedure relates 
a higher order moment to products of lower order moments 


and thereby effectively establishes an extra equation. 


ml these closure methods are used in fluid turbulence 
theory and so all might conceivably be acceptable candidates 
for closing a moment hierarchy in general relativity as well. 
if it is desired to describe cosmological models with 
a "small" degree of anisotropy and inhomogeneity, then it 
would seem that the appropriate truncation would be number 
two, the higher-order-moment-discard approximation. This 
emolce is reasonable because it would be expected that a 
inmement Of order ntl would be negligible compared to a moment 
of order n. Such weakly perturbed models are considered, 
for example, by Sachs and Wolfe [24] as previously stated. 
However, not all cosmological models are adequately 
iteseribed by small perturbations. In a recent paper by 
Woehler [34], a model is presented which assumes quantum 
Mmetuations of the metric tensor to be as large as the 
average metric tensor itself. Certainly the higher-order- 
moment-discard approximation is not adequate to apply to this 
model because all moments of all orders are important, espe- 
merely near the singularity. For this type of model, the 
memest Or third truncation alternative might be useful. 
Therefore, an attempt to construct a general perturba- 
tion technique based on a moment hierarchy scheme which is 
mele to treat fluctuations of all shapes and sizes is dis- 


eouraged by the fact that the method of truncation depends 


Lt 





Ween the type of fluctuation that is involved. ~iIn fact, 
=cn the solutions of the moment hierarchy could, in general; 
Mesquite different owing to the fact that different moment 
e@guations could result from different methods of truncation. 

5o it seems that a reasonable way to proceed is to 
eensider a model with a certain "kind" of perturbation 
field and then form and truncate its moment hierarchy by a 
technique that is compatible with the type of fluctuations 
feeumed. Thus, a two fold purpose would be achieved. First, 
a model would be constructed which describes a special non- 
uniformity about which conclusions can be drawn that may be 
Ceevalue to cosmology. Second, and perhaps more signifi- 
Cant , would be the “demonstravion that’ this procedure does 
Meum and 1s therefore relevant to perturbation problems 
eueeceneral relativity. 

Maeorder vo best serve the process of demonstrating 
Mmoemreeleyancy of the moment hierarchy technique to general 
metavavity, a Simple nonuniform model is used. It is Ssimpae 
maecme sense that only one function is required to describe 
meeeperturbation field rather than ten which are needed to 
meeat a general perturbation DrObl ena Nis. nmandomm une tion 
is also assumed to contain random perturbations which are 
small compared to the average field. Thus, the previously 
discussed higher-order-moment-discard approximation will 
temused for the truncation of the memene hierarchy which 
will be constructed in a way analogous to fluid turbulence 


Maeory. Even for this relatively simple model, the 
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Mavnemacical manipulations are quite extensive. It will be 
shown that this particular model is relevant to the study 

ero Jone wavelength perturbations and deformations of the 

Bsege Which approacnes the leneth characteristic of the 
PMemveroon ye sell. |) Thus. the model. Chough highly réestvicui ve. 


Mes a Special relevancy and direct applicability to 


Cosmology. 


aE, 





IV. METRIC TENSOR 


fee OFECTAL FORM FOR THE MODEL 
The process of building the model begins with the formu- 
tation of the metric tensor eres This can be accomplished 


byeeaving first the equation for the square of the line 


H 


element ds relative to a given coordinate system x and then 
computing Bee byeuhe 1elavion 
ae ak id 
ds“ = g dx"dx” . (4-1) 


mmr s end, consider first the line element for uniform 
Goomological models. Jt is called the Robertson-Walker line 


element and in spherical coordinates is given by [1] 





ds? = (ax®) - eb (XP) 5,2 (4-2) 
where 
do? = dr? + r?do? + r?sin?@d¢” (4-3) 
and Bea 
Bic 5) = ae ak Crs 
ry b + i — 
ro 


the parameter eh is an arbitrary normalization constant 
mipoecaamensions of length for the radial coordinate r. 

k is the curvature characteristic which may take only the 
memes —-1, 0, and +1. The time coordinate is x° and the 
meee Spatial coordinates are r, @, and ¢. The function 


RQ (x?) is called the scale length. 
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iober ti | points out that the spatial part of the fine 
element of Eq. (4-2) describes a three dimensional geometry 
feeety 1S homogeneous and isotropic with a uniform radius of 
eumrvavure Ra: Associated with the three values of the 
Puavature characteristic k are three intrinsically differ- 
Emeexinds of uniform geometries. The first, for which k 
memequal to zero, is the familiar Euclidean geometry. This 
eeometry is based upon the postulate that through any point 
mmmewoe Of a straight line, one, and only one parallel line 
Cam be drawn. The case where k is equal to +1 corresponds 
wo, vhne uniform geometry first discussed by Riemann which 
io) based on the postulate that any two straight lines ina 
plane always intersect each other. Finally, the case where 
k is equal to -1l corresponds to the noneuclidean geometry 
of Bolyai-Lobachevaki. This uniform geometry is based 
meermeune postulate that through a point outside of a straight 
mee, an unlimited number of straight lines can be drawn 
tem GO not intersect the given line and which lie between 
im@mcuralght lines that do intersect the given line at 
mmednaty [26]. 

Now, if the function G(x°,r) were replaced by an arbi- 
trary function M(x°,r,6,¢), the resulting line element would 
Meomlonicer be necessarily uniform. In fact, a positive 


function Ry (xr s6 5) could be defined such that 


Ry (%95P5659) 


oM(x°.r,65o) _ — 
2 r 
r, ; cE 7 : 


>) 





2 
an 


eal 





Was, at every point ae a tangent Robertson-Walker geometry 
@emlad be constructed with a radius of curvature Ra equal to 
mae function cu Nemec wat hewrune.: Lon Ry could be viewed 
jm this sense as a "local radius of curvature" or "local 
scale length." 

DoS tavedeprevloucgly One “Ol The Object lyes ol  Unas 
mmo serclation is to demonstrate the applicability of a prin- 
etple of fluid turbulence to perturbation problems of general 
memativity by means of a simple model. Clearly, a model 
whose geometry is described by the line element of Eq. (4-2) 
with G replaced by M is simple because the entire nonuni- 
formity is described by the single function M (in general, 
ten independent functions are needed to describe a nonuni- 
form metric B yx") because Bay is symmetric with respect 
jemu,Vv). Thus, such a model is a convenient choice for 


such a demonstration. It is therefore assumed that. the 


Tine element cf the model is given by 


0 
er ee os Nee. (4-6) 


Because this model is to be described by statistical 
meamoas, the value of the function M at any given point 
Inwee DOSSCSS a certain amount of randomness. That is, it 
must be viewed as a stochastic process or random field. 
momec, Mis a function not only of the four Space-Time 
coordinates (x’,r,8,¢) but also of an ensemble parameter & 
fimel: labels each of those functions M that are possible 


for a given experiment. Hence, Eq. (4-6) describes not just 


ee 





one geometry, but an entire family of geometries each of 
which must obey the Einstein field equations. 

ie MOtlLoOn of sano experiment in the study of stochastic 
processes is well defined. See, for example, Papoulis [18]. 
Mewever, for the purposes here, there is an interesting 
connection between an experiment (as used in statistical 
theory) and a modern concept called superspace. Wheeler 
[5] defines superspace as the manifold each of whose "points" 
meen esSelits one three-seometry. heaven: Submanitelid ol 
Bier space is therefore a history of the geometry undergone 
Mameelavseical space. Hence, the term experiment as used in 
pees paper would indicate an observation of the "state 
vector" of a possible three-geometry in superspace at a 
“—eemeuime. Cervainly, the set of possible state vectors 
for a given experiment is severely limited by Eq. (4-6). 

i’ the tensor form of the line element given by Ea. 


(4-1) is compared with the line element given by Eq. (4-6) 


with 
x0 = x0 
pate) =r 
X- = 64 
ee ae 


iets Casily found that the nonvanishing elements of the 
ieagete tensor for the model must be 
Eo = 1 


Gy =-e (4-8) 


GOntc:. 


aS 





f M 
G5, = —- r*sin*6e (4-8) 


B. DECOMPOSITION OF THE METRIC TENSOR INTO ITS AVERAGE 

AND FLUCTUATING COMPONENTS. 

In fluid turbulence theory, each of the random field 
variables (pressure and velocity) of the governing Navier- 
stokes equation is decomposed into its average and fluctuat- 
ing components as described in Chapter II. This same 
g@ecomposition may be performed upon the random metric tensor 
ety which is the field variable of general relativity. 


Waat 1s, en may be written as 


eae eG FF Sato ea) 


where en. is the average of a and Ks iS Sa Ue Lie wana 
eomponent . TorUseaveravce as UNdersvoode Cove aupensemu le 
arewage,. That 18, to take an average at each point xh, 
the Value of the random field from each member of the 
ensemble of possible states of the system which occurs 
Merl a certain probability are summed and divided by the 
memoer of ensemble members. AS this number tends towards 
infinity, this sum tends towards the average value. (The 
concept of average is discussed in more detail in Appendix B.) 
As seen in Eq. (4-8), the randomness of Bay is contained 
Mmemvice SLOchastic process ae. This process can also be 


decomposed at each point in space-time into two functions 


eechn that 
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M M 
e = <e > + Pv Heo) 


Nie 
where <e > is the average of ee and Py LSnenewr I Uelucae font on 
M 
about its average. By Eq. (4-8), the nonvanishing com- 


Bements of <g> and f must be 
uv LV 


<Bo > = 1 
<n) > = = eons 
<Bo > = - r2<el, 
2 Dee eee: 

Seo) eo Stns eae 
fa - - Pv 

ee 4 
f 59 ite Dv 
f 33 =a r*sin’of Cie dai) 


M 


MOeOCNet MSWOCNa Sue. Drocess t(x") Cai be Cet iMed- 25 sene 


moewo Of the fluctuation of a YO VtVekaverate. “Wnag te oss., 


ie 


eS : ‘ CHa) 
San 





ea i Pome xXDreSscea “AM Tero Ol Tart ls NOnVvanichine levenie 


become 
fo. = - eee 
2e4n2 M \ 
ewe cess in 65e" Sq" o (4-13) 


Because the average of f,, must be zero according to Eq. (4-10), 


M 
merlollows that 


<a> = 0 CS 


ao 





Also by Eq. (4-10), it follows that 
a ea (4-15) 


Beeause a and therefore <el'> MUSt De poSsivive, (Certain 
models with small perturbations would be ones for which the 
ratio of the nonvanishing components of ae Ze Eom as 
small; that is, for these models, the magnitude of +t would 
Eeamvery Small compared to one. As stated in Chapter 111, 
Weeers tnis kind of perturbation that will be considered in 
Paes work. 
ie RKequirement for the Average Component of the Metric 
rensor 
DUppOsSe Av iS-assumed that the average metric tensor 
“Ey? 18 Robertson-Walker. Then, the average geometry would 
be homogeneous and isotropic everywhere. This assumption, 
Moen is used extensively in the literature [7, 14, 24], 
aS convenient because of its relatively simple mathematical 
mom, Furthermore, in the limit of zero fluctuations, the 
model becomes determinate and the geometry is automatically 
Mm t orm. 
Haus, the nonvanishing components of the metric tensor 


“Buy? HOmme Ne enodel ares ¢ ily en iby 


<2o - = J 

Je aG(x" sr) 

ee SoG a) 

one o= re sintoed(% oP) (4-16) 
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Because the average metric is assumed to be everywhere homo- 
geneous and isotropic and is given by the above equations, 
mr tollows that ee = ae BY Virgue or Eqs wy (4=4) Sand 


(4-5), it can be easily seen that 


<Ryp (x° 47 ,056)> = Ry (x?) (4-17) 


where Ry (x? 575056) is the already discussed local scale 
temeth of the random metric and Ra is the universal scale 
Memeth of the average metric. If AR*(x°,r,0,¢) is the 


2 
fluctuation of Ry from Ra such that 


Ruan at AR oe (4-18) 


wmen by Eq. (4-12), it follows that 





(Ss ee (4-19) 


meus, t can be interpreted as the ratio of the fluctuation 
of the local radius of curvature squared to the average 
ieerus Of curvature squared. 
Eee ooouireneny tor the Kiiceua sine Component of ~vhe 

Mere lew Lonsor 

ier lua euurbulence theory, avcenvenitene-reduare— 
emo concerning the statistics of the perturbation field is 
Mee requirement of stationarity [13]. If this statistical 
@iletity is borrowed and applied to the present model, some 
simplifications can be expected by assuming the statistics 
of the turbulent field tT to be stationary with respect to 


a change of coordinates on a spatial hypersurface of constant 


a 





x®. Clearly, the statistics may not be assumed to be 


stationary with respect to the time coordinate x° because 
this would imply a model universe that ie not ina state 

of dynamic evolution (i.e. expanding). In Appendix C, the 
implication of this concept is discussed and the important 


formula given below is derived. 


9 .<h{t(x") ]> =O (12205 


where j = l, ae Syand W(t) is any function of Tierstochasrie 
Wmeecess t. Clearly this is a very useful formula for by 10 


-many bothersome terms are eliminated. 
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V. EINSTEIN FIELD EQUATIONS FOR THE FLUCTUATING METRIC 


A. DESCRIPTION OF THE FIELD EQUATIONS IN GENERAL 

Just as the Navier-Stokes equation relates the velocity 
ite hd vo the pressure field in a fluid, so do the Einstein 
field equations relate the geometry of space to the energy 
content of that space [1]. Being a tensor theory, the 
field equations are expressible in a form independent of any 
eeordinate system. The geometrical variable of the field 
equations has already been identified as the metric tensor 
ay The energy content variable is called the energy 
momentum tensor a which is developed and discussed later. 
ieeeshould be stated, however, that the field equavions for 
Pies model itself are not expressed in covariant form because 
a special coordinate system given by Eq. (4-7) has already 
been imposed. 

Blarting with the metric tensor eae Che [lela -equavions 
can be constructed as follows. (This procedure is given by 
Adler [1] and is used in the subsequent sections to build 
mame field equations specific to the metric tensor of the 
proposed model.) 


Horst. the metric tensor ae iS SUS bIGited AnGo ne 


meorlowing expression. 


dt 
= = + - -1l 
[uv 8] = 5 (a8, + 88, - 968, (oe) 
where al = anole} Pirosst 1S thevChrictol fel “symbot ol (une 
ox 
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first kind and is a set of 64 elements which are symmetric 
ieevMerespecy to their vt, Vv indices. 

Next, 42-quantity ie > LnervChristorrel Isvnibolsor tne 
second kind, which also is a set of 64 elements, is calcu- 


fated by 


Cm cas 
ee gE Btiy as] (5-2) 


where the identical indices indicate a summation over all 


an Bee, is also symmetric with respect to its lower indices. 
Mee object ae secalled- the second ‘rank Contravariant” Mepcre 


meme or. 1s related Lo the metric tensor Bie by 


ap O 
= 6 a 
E Eig u (5-3) 
Next, the Ricci tensor Roy iS. Shorined soy 
2 Oo p GA pP — p&% PP 7 
Ruy ore ic ae te it oie Tae We : (5-4) 


impeas a set of 16 elements. The Ricci tensor is then used 


memealculate the mixed Ricci tensor R* given by 
ene oe | (5-5) 


Next. the Ricci scalar R is formed by the contraction 


eeone mixed Ricci tensor R given by 
Pech ae (5-6) 


Finally, the mixed Ricci tensor of Eq. (5-5) and the 
maeei scalar of Eq. (5-6) are used to form the mixed Einstein 


tensor Gh 
Vv 
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ee 
Qo sR. -sorR . (5-7) 


The Einstein field equations are then formed by equating 


the mixed Einstein tensor oo to another tensor eT 


called 
v 
mie mixed Energy momentum tensor, an object which is related 
co iy aoGoloediseussed later. 
BS. the mixéd form of the Einstein field equations 


as they are used in this work are 


Tr = Gry (5-8) 
where 
ees g i 


EeeeCONSTRUCTION OF THE CHRISTOFFEL SYMBOLS OF THE SECOND 

KIND FOR THE MODEL 

mie task of constructing the field eaquations for the 
Bmoposed model itself may now begin as outlined in the 
Pee ecding section. From the field equations, the moment 
hierarchy will be determined. 

Pe iipOrrancoob ject. That. 1s immediately needed to -1form 
Wmemeooristotfel symbols is the second rank contravariant 


ifemri1c tensor gh V Tis col, ec &-can spe yoompu ced Tron has. 





fl=6) and (5-3). The nonvanishing elements of this tensor 
are 
one = ] 
pos - aM 
ans eu (5-9) 
E me Conus 


Syl 





=o Comes 


eo a Ceo) 
r*sin79 


The next step is to substitute the metric tensor of Eq. (4-8) 
into Eq. (5-1) which yields [uv,8], the Christoffel symbol 
emcne first kind. Because the second rank contravariant 


metric tensor gl Ls Giagonal. the calculataom of re the 


y 9 
Christoffel symbol of the second kind, by means of Eq. (5-2) 
is simplified. 

tee TeCatlculation of ee is still somewhat lengthy but is 
entirely straightforward. The nonvanishing components of 


i Eun Out LO DE 
LV 


Thi s = es om 

tS = ren dol 

Tae = = Benzo seme 

te = iv = - dM 

th = - aM 

te - a ol 

i a WE 7 = Jaye 

Te =e = r2(3 + =) 

Ty = - = r?sin?e(3M + <) 

Te = Ta = 5 90M (5-10) 


CONT. 
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r2 _ 11 COMi. 
a ae (5-10) 
r 
2a) Meee 2 al es 
Vote ato. 5 CaM - =D 
2 
Tn = - 34M 
2 2 
Pp = [43 = - a 
2 
T33 = = 5 sin*9(2 cote + 39M) 
3 3 al 
Y'o3 = Le rs D 0 oM 
ee 
Pes ae G 
3 3 il 2 
3) = 13 = 5 (3M a = 
3 at iL 
lr: = 0M 
° e sin29 ° 
3 = 2 = all 
3 1 
133 ai a oa 


mee CONSTRUCTION OF THE RICCI TENSOR FOR THE MODEL 

ime NeXt Suep icsbLhe eCalculavion Of the Nicel Tensor 
iy from Eqs. (5-4) and (5-10). However, it is anticipated 
meat Only the diagonal elements of the Ricci tensor are 
needed. (This will be justified in the subsequent section.) 


Meese diagonal elements are 


2 
Roy = > agM + @ (a,M) 
= 
Ry = 0M + = d_M + : : oa ea . 
a ee mesa oO 
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; i: eC aM) + = ; (om) 
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ee 


ae it 


LLCO Ue F 2 
> do M + 3- —— (9 4M) 


2 
r e-sin-o 


ee 

ai 
QW 
= 
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2M 1 
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r29°M + 3°M 
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=: 
Ieee = cee 6 


2 
d(M + 
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sin’6é 


2 ll 2 
34M + yr (3M) + 


iL 1 


alavede, 








NO] Re 


2 
(9 gM) 


le 
hope 


2 
Coug dM + 4 rol - z r2e(3 4M) 


2 2 
R33 = - = r*sin?6e" 35M 4S eine Meets an 


i 
2 2 


2 
+ 97M - é r?sin?6e"(2,M) + sin?6 coté 3.M 


> 6 


+2 wos oo Mm)? +i Sintecs M) 
y r q 6 


+ 
No] 


r sin?9a_M (5e40)5) 


WeeeecONSTRUCTION OF THE DIAGONAL ELEMENTS OF THE MIXED 
BINSTEIN TENSOR FOR THE MODEL 
iy order to construct the mixed Einstein tensor, the 
two objects RA and R must be found. 
dine mixed Ricci tensor R is easily calculated from 


d 


Homee(5>-5). Since gr tS vd lae0nal Tensor. kh . can be 


written as 


re) = 
R" = 6 weer (5-11) 


where the parenthetical indices indicate no summation with 
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respect to i}. Equation (5-11) justifies the assertion that 
only the diagonal components of ey are needed if only the 
eereaeonal elements of ae are Uo De calculated sbecalcewne 
me —Ciagzonal mixing occurs by virtue of the diagonality 
of or | 

The Ricci scalar R is easily calculated from the diago- 
nal elements of R* by Eq. (5-6). Finally, the components 
Ore the mixed Einstein tensor a piven by bg. = (o=/) (burn 


out to be 


ee 2 See ee Fae 
Go fF (aoM)* + e Lan + — 2gr 
2 
#—+— 9M + GT > (2,M) +f aM 
r2sin2g ° Tat 
2 zZ 
pace aM + 7p —— (aM) + i (3M) J 
r G-ein-6 
1 On 2 5 < —lM lila, 2) 
Gh, = - BoM - pF (BoM) + OTL 5 GM 
2 
+5 a t Sam + 5 SOke 9 om 
r*sin’6 r 
il 2 
+ T (3M) ] 
2 
C= - oe i (a.m) +e "[ 5 atm + 5 z 93M 
- Sane 
aL : lal r ices | 
+ ees (99M) + ee 3M + ; IgM J 
r r 
oes 2 3 2 Mp 1 42 ee 
Cea o“M - (dM) te [ 5 o<M + 5 2 34M 
(5-12) 
EOnes 
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ae 1 i ( .? COn te. 
7 _M ———— (9. M : (5-12) 
: y rm sin-6 ? 


+ 
NF 
S[k 


The off-diagonal components are not given because boey 
mee TOL needed. However, calculations similar to the ones 
required to obtain Eq. (5-12) reveal that the off-diagonal 
womponents of G are in general non zero and contain mixed 
partial derivatives of second order and second degree prod- 
Mets of first partial derivatives with respect to all four 
Seordinates. 

ite omould be noted that if Mis @-special function of 


maiy the coordinates x°® and r, namely 
M= G(x? ey) (5h) 


where G is given by Eq. (4-4), then the line element is 

identical to the line element of uniform cosmological models 

flJ. In this case, as should be expected, Eq. (5-12) is in 
r 


feu identical to the diagonal elements of G . for these 


eocmologies. 


B. CONSTRUCTION OF THE DYNAMIC EQUATIONS FOR THE MODEL VIA 
THE EINSTEIN TENSOR 
It should be noted that only two "dynamic equations" 
One oe are used to describe uniform cosmological models[1]. 
But each expression of Eq. (5-12) represents one side of 
feua “dynamic equations." If the proposed model of this 
Eumay iS to pass to a uniform cosmological model for arbi- 


Mearily small perturbations, then the number of equations 
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which should be treated as dynamical equations ought to be 
mao rather than four. 

In order to see how to construct two equations from 
equation set (5-12), it should be observed that G’, Contains 
second partial derivatives with respect to all four coordi- 
fMelces X°’,r,0@, and ¢. But each of the remaining objects Cate 
ee and G?, CCneain Second parcial derivatives: Witaeweopee s 
memonly two of the three possible spatial coordinates. 

This assymetry can be removed if Ce 1s taken "as one sidevor 


one of the desired two dynamical equations" and one third 


the sum ae, such that 
+ G° (522) 


is taken as the other. These two equations contain dynam- 
feat information concerning or and contain second partial 
@erivatives with respect to all coordinates. They are also 
found to degenerate into the usual two cosmological equa- 
tions when G is replaced by M [1]. 

Pwevircue OL Ha. (5-0), the two dynamic equations for 


Bae model are 


ae ahaa” (5215) 


and 


: qd 
: = (5-16) 


WIA 


From Eq. (5-12), a, is given by 


a a2M (5-17) 


e 2 _ 
ql. = - 392M - 2 (9M) +e “Laz 4 
J : : r cont. 
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7 1 2 ie 2 COM 
oe) at (0. M) + - — (3,.M) (Sato 
re sin-6 > age ae 
7 1 é 2 
ae SS Oe = (9 M) + eee a, MJ. 
r? sin?6 ? x ae 


For convenience, G can be expressed in terms of a", 


ao, = 6%, + 392M + 3 (90m) (5-18) 
KF. CONSTRUCTION OF THE ENERGY-MOMENTUM TENSOR FOR THE MODEL 
me proceed further, it is necessary to express the 
Prerey-momentum tensor in a form that is specifie to the 
teete used for this model. It is to be understood that 
the mixed energy-momentum tensor ee ts also -d stochastic 
mpeocess' because it depends, via the Einstein field equations, 
men the metric tensor which is a stochastic process. The 
Be@eral form of the enero to nonen eon Censor -1s. Known) Go 
depend upon three quantities, D the proper mass density, 
P the proper pressure, and u the velocity of the medium. 
These quantities D and P are termed "proper" because they 
mre referred to the proper frame which is that frame that is 
instantaneously moving with the medium. From any other 
ivame, the well known effects upon the mass and volume due 
memeehne relative velocity u must be taken into account. 

The representation of the general energy-momentum tensor 
meen respect to the metric that is used in this model is 
eemivyead in Appendix D. All of the elements of this tensor 
are given by Eq. (D-19). The two related objects which are 


important here are ue and te 5 
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Te tee = (Ratt) (5-19) 
Cc 


where 





D = proper density a 





l-u 

P = proper pressure P x — (5-20) 
ere 

UFeetocal velocity ol tie mavver. 12h t Ss 201 2c. 


beeause of the recent discovery of a possible aK black- 
body background radiation [6, 23], a relevant energy-momentum 
eemcor 1S one which describes disordered electromagnetic 
Maatation of this type. However, this radiation is believed 
to be no longer important to the dynamics of the universe 
(because now the matter density far exceeds the radiation 
density) although it is believed to have been in the more 
Gmebant past. The energy-momentum tensor which is relevant 
iemene current stage of expansion is one which describes the 
feamcon motion Of matter in the form of galaxies and clusters 
Semealaxies. Ea. (D-4) of Appendix D gives for these two 


Sars c S 


P = ac*D (S228) 


O if universe is matter dominated and 


where C 


a 1/3 if universe is radiation dominated. 
Hence, by Ea. (5-20), 


Pa nac ob). (5-22) 


BW 





Vi. MOMENT EQUATION HIERARCHY 


A. AVERAGE OF THE DYNAMIC EQUATIONS 

Equations (5-15) and (5-16) are the "equations of motion" 
for the assumed metric tensor of Eq. (4-8) and the energy- 
momentum CemMSOr “appropriaves Lor that. meucic. These equations 
should be viewed as the analogue of the Navier-Stokes equa- 
tion for V and P given by Eq. (A-1) because, as yet, the 
decomposed metric tensor of Eq. (4-9) has not been inserted 
mee them. 7 

memever, the substitution of the decomposed metric 
can be accomplished rather easily. sieoure xk inspection iO 
a, and Gee reveals that they are made up exclusively of 
metaeial derivatives of the stochastic process M and a factor 
a. meron Eqs. (4-10) and (4-12), it is easily found that 


= = Baal + T) (o=a5 


where 4 ie associated with the average metric and © with 
Pyemeructcuating part of the metric. 

Equation (6-1) can be used to separate a, and G’ TAGe 
two parts, one which contains only the nonrandom function 
emema the other which contains the random function T in 
@eonemnacions of Gand t. The part which is made up only of 
G terms must be , and Coe based on the average metric 
(because this part cannot depend on T) while the remaining 
part must be the source terms (because these terms must 


depend on T). The resulting expressions for en and q? . 


° J 
40 





are (using Eq. (4-4) which is the definition of e& and with 


me equal to Ro andaliGse: “Gh=i15 andr (5—10)))) 
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j " a 0 6R!' Soe 2 ona? 
ak = G = G - = - + = 5 
6R! Ot 
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it these equations are averaged, the result is the 
agesired dynamic equations for the average variables aoe 
aja) <T° >. 

These averaged equations can be immediately simplified 
Peyeusing the assumption of spatial stationarity that is 
eerseussed in Ch. IV. Thus, all averages of the form 


8, <h(t)> vanish by virtue of Eq. (4-20). Therefore, in 
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aru LCULAr:, 


| ; 
<as fy]? = 0. 


From this follow the identities 
Z 
2a at) at 
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and 


mus, with Eqs. (6-5) and (6-6), the average of Eqs. 


mracs (6-3) are 
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(6-8) 
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7 R? {)> Shane «ene, (6-8) 


If Eqs. (5-19) which give expressions Por ee and oe 
are used, Eqs. (6-7) and (6-8) become (with the aid of 


Ea. (5-22)) 
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These two equations are the desired dynamic equations. 
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Bee CONSTRUCTION OF THE SPECIAL MOMENT EQUATIONS 

As stated in Ch. II , it is necessary to determine the 
functional form of the extra source terms of Eqs. (6-9) 
and (6-10) in order to solve for the average quantities R, 
se > ; and <TY >. 

PMicwogctemdea 15/0 Construct additicenal equations 
mmeohe Unknown correlations. This is done by first multi-= 
plying the unaveraged Eqs. (6-2) and (6-3) or combinations 
of these by some "special functions" which involve powers 
Or derivatives of T and then averaging them. As might be 
expected, the nonlinearities of the unaveraged equations 
cause more unknown correlations to be formed which requires 
the formation of more equations. Although this process can 
Memmearried out ad infinitum, a special truncation of this 
meee sS 15 required which imposes an additional restriction 
upon the model. The truncation that will be used is the 
higher-order-moment discard approximation that is discussed 
mmeem. Lil in connection with the fluid turbulence hierarchy. 
moo means that the higher order moment that results from 
m@iewmonlinearity is discarded because of smallness. Hence, 
pmieeecne first "generation" of the moment hierarchy is 
meeced and it turns out that three special equations are 
sufficient (see Ch. VII). 

mine two special combinations of the unaveraged dynamic 


Eqs. (6-2) and (6-3) which are used as the basis for develop- 


imeechese three equations are 
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iemidus (6-11) 1s. mulgiplicd byvthe “Special function 
itgeecac, averaged, the Tirst moment equation is formed. The 
resulting equation can be simplified by again using Eq. 


(4-20) which implies that 
as<h(t)> = 0 


imeeect) is taken to be In(itt), then the following identity 


results 
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maus, the simplified first moment equation is 


2 
; (un 
keen Ia eet TY, )(1+t)> ss - eo") <| 


> 
Cita 


2 Z 
(on) (3,T) 
aes < ee eee > of Do eee < ae > 


+ : 5 
re (1+t) r?sin76é (1+T) 


SES uy 1 seeteo ae) (6-14) 


SO iC Oe 





If Eq. (6-12) is multiplied by 1 and then 9,71, the 


second and third moment equations are produced. 


2 
Corea (9 1) (rt) 
ech. = pd te a | + 2 |}. 


(9.7) (Tt) 
_ oR ( |; (6-15) 


a (1+t) 





cas eee) : 2,0) | 
K<(T ee + 3 )68ot)? = ac "C2 > i: a oe 


ie ae 
6R! : l rele: (6=— ioe 


oe Gees 


Again, using equation set (5-19), these three equations 


become 2 
(3 1) 
-G | r 


eet a eo ee) = = e =e 
ee) 








2 2] 
(3.7) (9.7) 


Peers r?sin?9 {(l+t)” 
6k 6R 6R! 
Sea Oe eee essa (6-17) 
R2 R?2 R 


(921) (1) 
K<[D(1+3a) +. Die (lea) err _ 3< ios > 


2 

(Oey i) Cont). 

+ SS: {| os oR (| a : = Q (6=15) 
7: (1+7) ‘ 
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3 2 
Chi) omar 
3 0 OR! | 0 | _ 
Pes. a > ane awe > = Q. (6-19) 
5 ae R (1+t) 
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VII. TRUNCATION OF THE MOMENT EQUATION HIERARCHY 


In both Eqs. (6-9) and (6-10) it can be seen that the 
source terms contain (l+t)” in their denominators. Because 
(1+1)” can be expanded in a power series in 1T, the source 
terms can be expressed as an infinite sum of moments of 
Various orders in t and its derivatives. As mentioned 
before, an equation for each of these moments can be gener- 
ated by multiplying the unaveraged Eqs. (6-2) and (6-3) 

Dy an appropriate quantity and averaging. The nonlinearities 
Oi these equations create a dependence of the desired moment 
meeorehner order moments. Because this occurs for all moments, 
gemethnod is needed to truncate the system in such a way 

mie Only a finite number of moments and required equations 
mete hierarchy are left. In Ch. III three possible trun- 
cation schemes are presented. In this work, the higher- 
order-moment discard method is chosen by assuming that 

ijest of the probability mass of the random variable cane) 
is concentrated around t = 0, So that. the naeher order 
moments can be neglected. This obviously limits the appli- 
mpi lity of the theory to models of small nonuniformity 

for which |t|<<1 (see Ch. IV-B). 

However, care must be exercised in affecting this trun- 
meron. Lt is not correct to just replace (1tt)” Dy ‘cone ran 
each of the denominators of the above mentioned equations 


must because this is the first term of the expansion. This 


48 





can be seen by an example. Consider the typical term 








{os |> ‘ Thais term can be expanded AS. 1 OleLows 
ET 
d)T 
eae oe eae ek a) ater, C7) 
Since ees oe Toe 0, 
ore 
Saleen |: we chante) ole (yee) 


However, if (l+r) were falsely replaced by one, the follow- 
ing incorrect conclusion would result which disagrees with 
mae correct Eq. (7-2) 


erat 
< area > es “oe ; (7-3) 


mence, an expansion must be made for every term and that 
mowest order moment term which does not vanish by virtue 
of having:a zero average is assumed to be a "good" approxi- 
ie On CO the term itself. 

With this "weak fluctuation" assumption, Eas. (6-9) 


(6-10), (6-17), (6-18), and (6-19) become 


2 Be AN 2 3 an 
K<eD(ltau*)> = > a4 <t ee Heo) > 
os 2 2 2 
= 2 e G ON eet: = <(3,T) > + : <C947) > 
r r? r*sin76 
2 
1 
+ 3k <tT*> = - 3k _ 3K) (7-4) 
R? R* R* 
K ' ny 
~ <Du2(2a-1) - 3aD> - 2B ajcr?> - 5 at<t?> (7-5) 
3 Cone... 





3 5 =e 2 1 2 
+ i Sea) a Ts e }<(,.0) > + ae! > 
+ : < 
CUO Ek oe ee eS 
eosin > fe 
2 
\ ORM R? Conc. 


-- = - 5-2 (7-5) 


ae 2 2 2 
~ ze [cag > + - 1055) Pa —— at > 
r r°sin’6 


+ 2 K<D(1-u?) (1-30) (1+1) 


Z 
Se Get, Zone, (726) 


K<[D(1+30) + Du?(l-a)]ro- = ace 


; | 2 
- G26 > Su Cera )s = (7-7) 


0 


2 2 
3<(8,T) > + aRt <(3,T) > 


= K<[{D(1+3a) + Du’ (1-a)]a,T> (7-8) 


Lope 


im the limit of zero fluctuations, Eas. (7-4) and (7-5) 
become the usual two cosmological equations for p and R 
that are given by Adler [1]. Thus, it is natural to look 
Upon these equations as dynamic equations for o and R. 
fne source terms which appear in these equations contain 
the three functions <t?> , <(a,7) >, and (<(a.1) > 


1 g i ° 
+ ar <(,7) a cia <(3 47) Sie, Unvsy only Thre 
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equations are needed to determine these variables which 

justifies the statement made in Part B of the preceeding 
chapter that only three equations of the hierarchy need 

to be calculated. These three equations are Eqs. (7-6), 
(7-7), and (7-8). Although the above functions <t*>, 

2 2 ci 2 1 2 
o,t) >, and {< (9,7) > + aa <(9,T) 2 Ee peer <(84T) > 
Bepear in these equations, additional moments also occur 
aoreh relate the density D to t and Dtod,t. It is shown 
ma the succeeding chapters that these moments can be related 
Memthe existing moment <t*>. Thus, Eqs. (7-4) to (7-8) 
are soluble among themselves. Because Eqs. (7-6) to (7-8) 
mevyern the dependence of the source moments on R, they 
mee called growth equations to convey the idea that they 


determine the development or growth of a perturbation T with 


maecreasing R. 
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VIII. REVIEW OF THE RESTRICTIVE ASSUMPTIONS 


Berore spyocecding further, it might be heloful to review 
the restrictive assumptions that are used to obtain Eqs. 
(7-4) through (7-8). The first is the very restrictive 
mom Of the metric tensor. This form excludes many types of 
gravitational wave phenomena. All that is seen by an observ- 
er would be a local "expansion or "contraction” of Space 
Because, as stated in Ch. IV, it would always be possible 
Pemconstruct an expanding or contracting Robertson-Walker 
geometry that is "tangent" to the given space at any given 
point in space-time. xe Oxcdiusi1 ON -O. “Succi awWaves imo 
Seavainly exclude any hope of using this model to describe, 
fmma realistic way, the violent, multiply connected space 
iors «Wheeler envisions must typify the quantum fluctuations 
of space-time. 

MaemSSCONd andebhiard PESLrICcCeylOns are the Tres~vricuilons 
Concerning the Robertson-Walker form of the average metric 
meoe the spatial stationarity of the fluctuating part. These 
evo assumptions rely more on intuition than upon scientific 
imeor. They imply that there exists some sort of Spatial 
Pmt ormity in the statistics of the assumed fluctuating 
geometry. However, it is to be understood that a "new" 
cosmological principle is not being seriously proposed here 
by the use of these two assumptions. That these assumptions 


are used at all reflects a need to simplify the mathematics 
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of the moment equations in some reasonable way by appealing 
to some plausible symmetry arguments. Thus, these two re- 
Suractions should be viewed only as steps in this particular 
preoolem whose purpose is to demonstrate a technique. 
Miia. the LoOurth rescvriction -is the assumption 
meat the higher order moments are negligible owing to an 
assumed smallness in the anisotropy and inhomogeneity. This 
aeaoumption 1S a restriction upon the applicability of the 
ween, that is, in those solution domains for which the 
Perm vude of tT 1s comparable or larger than one, the use of 
the model is unwarrantable. If large amplitude fluctuations 
eles LO be considered, then perhaps a truncation based upon 
pmomenomenological or a cumulant-discard approximation would 
De useful. However, in the case of the latter, it might 
memassumed that the probability distribution of the random 
weresadle t is normal. Then, the higher order moments can 
be shown to be products of the lower order ones [18]. Such 
an assumption is made in fluid turbulence theory [3]. But, 
meine case of this model, T must be greater than or equal 
memminus one as given by Eq. (4-15). This creates difficulty 
With a normal distribution where the probability of t being 
timers equal to the probability of it being -T where T is an 


feeicrary positive number. 
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IX. SPECIAL UNIVERSE MODELS 


Erem thas point on it 1S convenient to consider two 
timaoes Ol models which correspond to two values of a. ‘The 
feese, tor which © is equal to one third, corresponds to a 
universe made out of electromagnetic radiation while the 
Second corresponds to a universe which consists of cold 
mMatbcer for which a is zero. (a is equal to 


22 ).) 


a See Bae 

It is not practical to treat these models simultaneously 
by leaving ao arbitrary. The cold dust model can be solved 
exactly owing to the fact that all matter moves freely along 
geodesics without being accelerated by any pressure gradi- 
ents (pressure LS, Zeu0 everywhere). But the radiation model 
G@ees involve pressure gradients which do not allow the 
elements of the radiation fluid to move freely in space [9]. 
iemetretore, these two types of models are sufficiently differ- 


Saeeirom one another to merit separate consideration. 


Pee RADIATION FILLED UNIVERSE MODEL 
For the radiation filled universe, Eqs. (7-4), (7-5), 


eo), (7-7), and (7-8) become 


2 = 2 
oa _ ote : Sgt = e fet > 
2 
+ <(3 as es < (ah) | (oat) 
ee 0 Pood o ? Cont. 
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where 


and 


2 
. cont. 
4, 3K Be AN os. Sie (9-1) 
R2 R2 Be 
K 3 R! 1 : 
es = > RO aps aot Pe + ; <(3 7) > 
_ 5 .-G 2 1 é 
75 ¢ {<(aen nie = <(3,T) > 
1 2 
+ ———- <(9 +t) >} + so <p> 
r?sin2e@ ? ne 
: kK ei! Ri 
ge (teens Se (9-2) 


. i ofan)" oo <9 t) > —— <(9 05] 


2 
_~-L5k TRG eae 





=S25 - 2 - (9-3) 
Ro R? R 
Piccitd>) = 2 eG ae = ve ieee 
2 
=i SG) > (9-4) 
2 t zZ 
9y<(8yt) > + EB cage) > = dxc(ajr)a> (9-5) 
2 
pee =e (9-6) 
Z 
dies) Plat oe =20 (9-7) 


ieeonOould be noted that P is equal to the average of fa 


femmer [1] points out that ee can be interpreted as the 


Veval mass density with respect to the given coordinate 
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eyocvem and therefore includes relativistic corrections to 
mee rproper mass density. 
1. Elimination of the Spatial Derivatives 

NOweconsider the terms which occur in Eos. (9-1) 
and (9-2) of the form <(at) >. It should be noted that 
eauly 9, and on are’ partial’ Gerivaltives with respeccrtosa 
length (in the time and r directions). The remaining par- 
tial derivatives 9 6 and o4 are partial derivatives with 
mespect to the coordinates 8 and $¢ which are not lengths 


tmoechemselves. An inspection of the line element of Eq. 


a Ele» _ Ge G/e 


(4-2) suggests that 3,, . 


l/r 3, and e- 
irs ind 84 all of which have dimensions of inverse lengths, 
play the same roles in the x ee Ol and @ directions respec- 
mmrely, Because disturbances of the metric propagate with 
the speed of light [11, 32], it would be expected that the 
Characteristic inverse lengths associated with these quanti- 


ites are comparable in magnitude. Thus, it seems reasonable 


imomexpect that the spatial terms 


2 -G 
e%<(a,t) >, —= 


re 


=C 2 : 2 
a T > 
e Pay, ree , ( ‘ ) 


Nh 


(r sin®@) 


oS 


could be expressed as mere repetitions of the temporal term 
<(3,t) > multiplied by a paseeucdent of order one. 

Such a fact would simplify Eqs. (9-1) and (9-2) 
because the terms involving the spatial pertial derivatives 


would be changed into the already existing terms involving 


the time derivative. 
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This COnsOLIGation of —the Cerms is performed in 
Appendix F. Equation (9-3) is used to eliminate the spatial 
partial derivative terms from Eqs. (9-1) and (9-2). An 
eppropriate renormalization is then performed to bring the 
mene hand sides of the resulting equations back to "canoni- 
cal form" as they appear in the right hand sides of Eqs. 
(9-1) and (9-2). Equation (9-4) is used to change the form 
of the "renormalized" equations into a more convenient form. 
mimes last step is not part of the consolidation itself 
eemce it just allows the resulting equations to be presented 
mee more desirable form. The transformed pair of Eqs. 
eo=wee and (9-2) which are void of spatial partial deriva- 


Tives are 


Z ! 
out z K<tTd> + iz eer) > 5 3 ~ ee 


2 
R! 
- 2K ..2, er om eis (9-8) 
Re Re Ke 
2 
eS oe = i e<td> - 5 Senta + « — oe 
2 
lee 2 k Ga ate Fea 
SRere COD Re ee ee © ee (9-9) 
Re R* R R? 


oe inclusion of the Dynamics for che Ractavronne 416.0 


Because Eq. (9-3) has been used to remove the un- 


known 


e 2 2 iI 2 
e : <(9 t) > + = <(39.t1) >» *+ —— «(3 Tt) > 
a r2 0 r*sin’é@ > 


from Eqs. (9-1) and (9-2), there exists now only four 
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equations in the system namely, Eqs. (9-4), (9-5), (9-8), 
and (9-9). But this system of equations is insoluble because 
miley represent four coupled nonlinear differential equations 
in six unknown functions of the time coordinate x’. These 
MigiOWN Iunctions are op, R, ean <(a.t) >, <td> and <(0,t)d>. 
However, there is a way to express <td> and <(0,1)d> 
memeeerms Of the other four unknown functions. It is not 
unreasonable to expect that there should exist some rela- 
tionship between the density fluctuation d(x") and geometri- 
cal fluctuation +(x") because in some sense the geometry 
memecuations "drive™ the density fluctuations. But it must 
memremembered that the radiation fluid possesSes internal 
meessure forces which push the fluid elements about as 
well as do work against any expansi0ns CY Cconlvracciions oF 
space which might occur with amplitude + [9]. Thus, the 
MmeoL.em Of determining this relationship between d and T is 
immeceneral very complicated. The correct way to do this 


mmeo use the equations of motion for the radiation fluid 


which are given in tensor form by [1] 


sae mn cea (9= No) 
where 


UV, & UV yg pV pho 4 pH mov 
mo ott oe Day 


inmcieed, the relationship between d and Tt is approximated 
in Appendix I via these equations and is expressed in terms 


ei the correlations <td> and <(0,1)d>. 
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An instructive yet simpler way to determine the rela- 
Mbenship between d and t is the following. Consider a small 
SB-volume element V at rest with respect to the coordinates 
Sano, ¢). Let M be the total relativistic mass content 
wmeunes 3-vOolume. The total relativistic mass density D 


jmeecnus given by 
_ i 
D 5 V . (9-11) 
Peeemiatl change AD of D can be in general a result of a small 


Pmanee AM of M and a small change AV of V. That is, 


pp = ai _ Mav 
V2 





(9-12) 


Nowsthe quantity we td ‘x vis known torbe a “scalar 
Mmevartant where g is the determinent of the matrix Buy 
[2, 32]. In particular, from special relativity, this quan- 
imuey is the product of the proper time interval and the 
proper volume. Because x° has already been identified as 


the time coordinate (see Ch. IV), the volume V must be 


=v ood = els? 
faecerms of the coordinates r, 9, $¢, 

V = Y-g drdédd . (9-14) 
The use of Eq. (4-8) gives 

v = e242 (ar) (rae) (r sineds). (9-15) 


0 Ue eo iae) (9-16) 
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so that V = V¥(1l4r)2/2 (9-17) 


where V= @iG/e 


(dr)(rd@)(r sinédd¢) 

ine change of mass AM can occur because of the work 
mome Dy the pressure against the change of volume and by 
ime lt lux of mass out of the volume. If this flux is tempo- 


rarily ignored (i.e. all mass moves along geodesics), then 


miecan be expressed as 


AM = - Pav (Oo =4,39) 
e* 


mimere P is the total relativistic pressure. The relation 


p= 402d (9-19) 


AM = - SN Sp. (9-20) 


Then Eq. (9-12) becomes, using Eq. (9-20) 


ee ga 3] 
AD = - 3D >; (9-21) 

But Eq. (9-17) with small becomes 
veVistr. (G=22) 


ijmey is replaced by 2 Vr, Vo bya Vv ab by: dy. and. Db iby pan 


Eq. (9-21), then this equation becomes 
ao | (9-23) 


Hence, the correlation function <td> and <(9,1)G> become 
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<td = =) 2o<T- > 
and 
<(a,t)d> = - pd,<t*> . (e207 


Mime ppendax 1, the effect of “mass leakage" from the volume 
is approximated by means of two variable parameters A and 8. 
it is shown that a reasonable form for the correlations 


ard> and vena) de are 
<td> = - 2)p<t*> 
SGC UG Pa Bp Oates 2 (9-25) 


The coefficients A and 6 may be thought of as representing 
a zeroth approximation of those functions of x°, which 
correct the zero leakage correlation Eqs. (9-24). These 
eeemracients, then, become open parameters which character- 
iena given model. Certainly, if mass leakage were nov 
iaoorbant, then A and 6 would both be one. 

The substitution of Eqs. (9-25) into the dynamic 


Eqs. (9-8) and (9-9) and into the growth Eqs. (9-4) and 


(9-5) gives 
2 1 
KO - = Micon po oat 2 SO, Bot os ue 
2 
1 
~~ ct? oy ato (9-26) 
Re R* Ro 
K ‘ ah de 2 
eae KR - TZ <(997) 2 vores > 
Z 
eg eg otra a ha ee ae (20) 
R? Re R ie 





J2<72%> eae eer ee 2 = 2 : 
cr R ne > 3 MEOSt. = <(Qt) > (9-28) 


2 WR 2 y 
det) an a <(9 1) ae ai e Boe pai sen (9-29) 


3. Solution of the Moment Equations which Govern Growth 

Sheek TNeulaGLoOns 

The four Eqs. (9-26) through (9-29) now contain the 
mour variables p, R, <t°>, and <(31) > anor Nessy Suem ia 
now closed. However, the solution of these four equations 
isenot easy because they are nonlinear and nonlinear prob- 
memos dO not in general yield easily to analysis. But it 
foe possible to take advantage of the requirement that the 
moments which appear in Eqs. (9-26) and (9-27) be small 
mempared to Kp. That is, the additional stress terms cause 
emomalil perturbation of the usual, fiuctuation free cosmol- 


Ogical solutions for p and R which are well known. These 


seeutvions, derived in Appendix E, are 





Ko = - = (9-30) 
R 
and 
2 a 
ee so k (9-31) 
R 


where a, is y. Constant Of 1nvecrat leon and. kis Lae -curvavure 
@eeracteristic and may take one of three values -l, 0, #1 
(see Ch. IV). 

Nonlinear problems can in many cases be Samp iat Led 
by the replacement of a variable with a very close approxi-— 


mation to that variable [30]. If such a replacement is 
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mace fOr R!' and p in kas. (9-28) and (9-29) with the fiue— 
fmmaci1On iree Solutions Eqs. (9-30) and (9-31) and if ail 
derivatives with respect to x® are transformed into deriva- 
Mmeves With respect to R by the chain rule for differentia— 
tion, these two equations become soluble between themselves. 
However, as can be seen from Eq. (9-31), the effect 
of k upon Re is not significant unless R is large enough 
poe chat a/R is of the same order of magnitude as k. It 
is well understood that that value of R for which the effect 
of k is important to the dynamics is the current value of 
heme For this reason, much research in observational cosmol-— 
Oey 1S devoted to the determination of k on the basis of 
data on current expansion characteristics [25]. But the 
eoolaicability of the radiation model is only valid up to 
wae time when the radiation density is equal to the matter 
density for which the associated R-value is much less than 
the current value of R by several orders of magnitude [8]. 
Mens, if the R-domain of the solutions <1*>(R) and 
(),T) > (R) 1S restricted towalues Of i SUIT UCtenaly Less 
than the current value of R so that Sa |p then the 
following substitutions for R' and op in Eas. (9-28) and 


@5o=29) are valid. 





3a. 
KO a 5 (9-32) 
R 
and 
2 a 
Ria eaeee. (9-33) 
R? 
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It should be noted that the neglect of k in cosmol- 
ogical problems is not new. k is neglected, for example, 
byeHoyle and others in their expression for R'* for exactly 
the same reason that it is neglected here [10]. 


Wae erow ch equations. become, USing the above appro, 


imations 

Roa <7 2s + RA_<12> 8ir<tTe> = 2R! (3 _ ( 4) 

dp i R aL a iv = A < o | - a=) 
Gree (9.1) 
—Coet ES 4 (G7 ate ee 
0 L 0 4g 2 
5 A; oa eee Oe ce = — 0.<T*>. (9-35) 
R a, R da. RS R 


These equations are solved in Appendix J for <1*> 


2 
ana <(a,T) >. The solutions are 


i D , 
= 1/3 p,it4 / 3 tenga 
<-> .= B23 + Ses + 13 R (9-36) 
Soe) > My 4 6 
z = J _ 
: ae UP 3k + Ty /3R” + Dy/3R (9-37) 

where 5 = —- 4 + 27%2(0\+B) 

j= -4- 2V2(X+8) 
and 1473 , Di 73 ; B13 are arbitrary constants of integration. 


4. Special Submodels 


ieee Da Scuos Lor 
In general, a given combination of the parameters 


and F specify a nonuniform model which 


m3 * 1/3 ° 1/3 
is difficult to treat analytically. Yet, there are special 


values of these parameters which designate models that are 
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relatively easy to treat and which are interes vine Pinese 
are the three models for which only one of the three arbi- 


merny paramevers are non zero. That is, 


z 
Gir) <t*> =". fF . Shoe = hy F -4 
WS. as, re aw (9-38) 
I Ce 
: <(0,T > 
(2) <t?> = 2d pi | Com (9-39) 
hy 
D oe 
aes < T > 
(3) <1t2> = a5" ae —— = ae . (caer 


rt 


petually, not much generality is lost by considering the 


above three models separately. The I and F 


ee Ve 
terms of the general Faqs. (9-36) and (9=37) sare functions 

of R to a power. In general, however, these powers are 

@uate different from one another (i.e. try A, B equal to 

one which are the zero leakage values for X} and 8). Because 
mechanges by many orders of magnitudes during expansions 
eaemconcractions, only one of the terms would dominate 
anyway (except for those "particular" values of R for which 
two or even three of the terms would be of the same order). 
Hence, most of the models with a given combination of I, 7/3 


and F 3 Harameters woulda cegzenerate into One of the 


ma 1/ 


pgeve three models for almost all R-values. ’ 


However, there is a difficulty with the ya 


model. Equation (9-25) says 


ae tee) 
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X must therefore be a parameter which measures the response 
of the radiation fluid to a given geometric perturbation. 

It seems unreasonable to suppose that the fluctuation d of 
the total mass density is positive for positive +. That is; 
it seems that only models that are physically reasonable 

mee Uhe ones for which i is positive. But, from Eq. (9-38), 
only the negative A values permit the obviously positive 
fume tions <(3,t) > and <1t*> to bé positive a PSm Ost Give 


byevirtue of Eq. (9-32)). Hence, the F -model does not 


17S 
seem to have much physical meaning and is therefore not 
considered in this work. (If some positive \ values would 
Bewacceptable then there would be grounds for exploring 


the content of the F -model. ) 


LAS 


The two remaining models, namely the I and 


nS 
Dy ,3-models , are very different from one another. The 
Ty 37model represents a model which becomes more nonuniform 


with increasing R while the D —-—model depicts a model for 


Ayes 
mmrem the nonuniformity decreases with increasing R. 

Because there is no intrinsic reason why positive A values 
maememnou acceptable, both of these models can be considered. 


(I is used to denote the arbitrary constant of an 


7 3 
Mi@emcasing type perturbation field of the metric for a 
maddation, @ equal 1/3 universe and Ds 73 is used to denote 
Pmeedecreasing type. The same nomenclature will be used 


to specify the increasing and decreasing type perturbations 


of the dust, a equal zero universe but with a zero subscript.) 
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Ut ew Ol ANverecst. to see Noverne Tilvervuaticn 
peresoses Which correspond to Doth the 1,73 and Dy s3~type 
Pervurbations act upon the average total mass density o0 
mae the average scale length R. This can be seen if the 
stress terms which appear with the average density o in 
Fas. (9-26) and (9-27) are expressed in a form which cor- 
meoeoonds to Che functional form of the T1473 and Ds 37 type 
fluctuations given by Eqs. (9-39) and (9-40). Equations 

k 


(9-26) and (9-27) become (with ae <t*> neglected to 


Ko<t*> for reasons already stated) 


(I, ;3-fluctuations ) 
Ri [60 + 248 + 6Y2(A+8) |] 


Ko + C(I, 73 72 
Z 
eh ee OS toa Oe Ne, (Ge) 
Re Be 
Ke i £12. + 88 - 2V20+8)) 
a. 72 
ae 
ee ee an eee (Gees 
R* R Re 
(D, ,2-fluctuations) 
i x De eee iy Oe) 
ae 6(D, ;4)R" [60 + - (A+B ) } 
2 
eee oe ee (9-43) 
R? R? 
K 3 [12. + 8B + 2v2(A +8) ] 
eee aes 72 
2 
Se ee en (9-44) 
R? R Re 





our a oD a 
= ee - Nea 
ee. fy eg OLD) 


i= - 4 + 2720448) and j = - 4 - 2V2T)+8) . (9-45) 


Ttysnould- ibe poOinted out that of (the three 
parameters C(A), »4, and B of Eq. (9-45) (where A represents 
Ty 3 or Dy 73)> only C, lends itself to physical measurement. 
Heemust not be forgotten that Ch is proportional to a para- 
meter which is the result of an integration and therefore 
is an arbitrary constant which must be fixed by some obser- 
weer On. But A} and gp are not the result of any integration. 
Bbermeoare Variable parameters which, in a sense, measure the 
Mmeenematical uncertainty between two pairs of correlation 
functions (see Eas. (9-25)). If it were known how to 
exactly determine these relations, then 8 and A would be 
known numbers and hence be of nonvariable status. As it 
is, they have to exist as variable parameters. (In the 
dust model considered later, the dynamics of the matter 
are determined exactly. Hence, the correlation functions 
can be found precisely and i and 8B do in fact, for that 
case, become known numbers.) 

The solution of the dynamic Eqs. (9-41) to 
fo=44) are given in the next section. 

DewcC tessa, eau lon 

The question naturally arises: What is the 
effect of a given combination of fluctuation parameters 
C(A), 4, and B upon p and R? It has already been pointed 


out in Section 3 of this chapter that the effect of the 
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fluctuation stresses upon the average density p and scale 
length Ris small. Thus, the actual solutions are small 
perturbations of the fluctuation free solutions and it 
therefore seems appropriate to compoare the solutions of 
Eqs. (9-41) through (9-44) for p and R with these known 
fee uuation free solutions. 

This can be done in the following manner: There 
mee Only four distinct ways that the fluctuation stresses 
mieeperturb the solutions for 9 and R from their fluctuation 
meee SOlutions. That is, at a given time a the density 
op and the scale length R can each be perturbed above or 
Below their fluctuation free values. Let each of these 
four type of perturbations be designated by the letters K, 
meee’, and N in the manner shown in Fig. 1 for the increasing 


ir -fluctuations and Fig. 2 for the decreasing Dy 3-tluctu- 


E/ 3 
mgmoms. in both cases, a given letter indicates whether 
ie po and R are greater or smaller than the fluctuation 
ieee SOlutions for 9 and R respectively. For example, a 
model for which o(x®) and R(x°) are both greater than their 
respective unperturbed values p*(x°) and R*® would be clas- 
Sitied as type M. 

The manner in which the solutions 9 and R are 
perturbed from their fluctuation free values which are 
specified by the above simple classification scheme is 
determined solely by the two parameters werand 6. Ihe mag 


nitude by which a perturbation of a given type occurs is 


determined by the parameter GUA) A -Phise Vane Bed etermine 
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the qualitative behavior of the model while C(A) determines 


ime Guantitative behavior. 


Parerder to determine the classimiueataomonr eo 
model which is characterized by the three parameters CCA), 
A, and 8 where A refers to either the increasing I) /37type 
Peeourbation or the decreasing D137 UYPE the solutions “lor 


peoma KR are needed. In Appendix L the solutions for Ko 
2 
and R' are given in terms of some special parameters Py 


Da» and My which are functions of only the three parameters 


C(A), A, and 6B. These solutions are 


m 
ae 3a,R- aan e (9-46) 


Rieroene ys tebe (9-47) 


The unperturbed solutions 9* and R* are those 


mere which Py and by are DOL ecual torgero, ‘trp. 16 east 


een that the direction of displacement of po from p* is 


eewermined solely by the sign of f The sign of the para- 


A° 


2 
ewer b Onetue Obhereiand. “Cells tywnetner hs -l5starwcer 


ee 
x 
Sumsialler than its unperturbed value (R'§ for all values 


Maer If, in the case of the 1 —Yype model, a adorvwerd 


is 


time integration of Eq. (9-47) for various values of - 
jmommade from R equal zero, a family of solutions R(x?) 

@eme be generated. The displacement of any one of these 
solutions from R* depends on the ratio Ri? yer)”, For example, 


femchnis ratio is greater than one for all R-values, then R 


is necessarily perturbed above R* because all solutions of 
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the family are "pinned" to the same initial value. Because 


the size of this ratio compared to one is determined by by 


meetOllows that Ta 3 also determines the direction of the 


ay ae 


feesplacement of R from R*. It can be shown that the sign 

on D173 also determines the displacement of the decreasing 
Ds 73 solution of R from R* if the integration is begun 

at a sufficiently large R and performed backwards in time. 
in@eappendix L, the signs of fh, and by (and hence the classi- 
fication letters) are determined for all values of \ and 8 
mor both the D5 /3 and lip Uyeoe a LuctUuatlions.. Whesresiite 
are graphically displayed in Figs. 3 and 4. 

It might be asked how a given combination of A 
mde affect the exponent of R in the fluctuation stress 
em. Lhe agependence of this exponent is presented in Fig. 
meeierra plot of the lines of constant j and i in A-8 space. 

This ‘concludes the basic discussion and «calcu— 
meeons for the radiation model. In the next section, a 
Similar treatment is made for the dust model. The discus- 
S10n and comparison of these two models are given in 


tmapcter A. 


Bee DUST FILLED UNIVERSE MODEL 

Some of the calculations and discussions in this section 
mae similar to ones of the previous section which deals 
mabth the radiation filled universe for which @ is equal to 
mmeethird. Here the model for which a is equal to zero is 


treated. It is assumed that the reader is already familiar 


fa 





with the previous section so that otherwise redundant expla- 
MmaciLons may be omitted. 
For the pressure free universe, Eqs. (7-4), (7-5), (7-6), 


(7-7), and (7-8) become 


’ 


2 
Kp - 2 on ae z <(9 7) > 


7 ft 8h eCapa) + a EC aes oa <(9 05] 


r 2 6 r*sin*é ? 
2 
ye cai eee ee (9-48) 
R2 R? R? 


R 4. 3 a 
5 ie > es a cia +7 Sitges) 


I 
Rofo 


6 p 


= 2 2 2 
~ ee Ne(ait) > + = <(a.t) > + ——— <(a,7) > 
= r? i sane oe 


| 1 
+ eee = . a ee soe (9-49) 
R? Re R R? 
2 2 2 
-fe- a5 t) > += <(at) > + . COST) ac 
r r2 r-sin6 ° 


2 
LK ue Lok 





5 5 
=-—- ={K9 - = K<Td> - == - = - = (9-50) 
e 2 R? ‘R? R 
! 2 
K<td> - = Gee - oon a,<t > = - Boones > (9-51) 
2 ! 2 
a<(8jt) > +e <(agt) > =F K<(9,T)A (9-52) 
where 
Oe > 
and 
gd=iD-o0 =. (9-53) 


re 





ft should be noted that for a equal to Zero, the pressure 
Homiaentvically zero. Hence, there are no pressure gracients 
poyeperabie CO accelerate the matter. This implies that 
Mawcer, when left alone, will follow geodesic lines. If 
mmese geodesics are taken to be geodesics of zero velocity 
then the matter comoves with the coordinate system (r,6,¢) 
and the above equations follow. From Eqs.(9-53) it is seen 
that p is the average of the proper density D because the 
relative velocity u is zero in which case D is equal to D 
(see Eq. (5-20)). 

Mone tiie Ubon Ol. one opaulal Derivavives 

In Eqs. (9-48) and (9-49), the terms which contain 

eeavlal partial derivatives may be eliminated explicitly in 
amanner analogous to the way they are removed from Eqs. 
@—-1) and (9-2). This computation is performed in Appendix 


foeeelne resulting two equations are 


Ps ! 
Keo z <0. 0) > + # ee + 2 K<TC > 
2 
! 
he z KU dp ae ei fo See (9-54) 
Re Re R? 
a 1 R! io 
ST Cote > + o RO ond = oe lp Ea a> 
2 
tf { 
~Flar-.-f-8 -i (9-55) 
Re Re R R? 


a inelusion of the Dynamics for the Matter 


In order to solve the system of equations, it is 


necessary to express <td> in terms of existing correlation 


cS 





mIeoLONS. This 18 also done in the PFeV ONS) see i1on) oop 

is complicated by the existence of nonzero velocities. 

For the pressure free model this calculation is much Casicrs 
Consider Eq. (9-12). Because the velocities and 

the pressure are zero, there is no mechanism to change the 

mass content of an expanding volume. Hence, AM must be 


zero and Eq. (9-12) becomes 


AD = oe (9-56) 


With Eq. (9~22) for small tT, Eq. (9-56) becomes 


(also with the usual definitions for the increments as in 


the previous section) 


oa (9-57) 


Ou 
i 
rojo 


lemee, <td> and <(?,1)d> must be 


<td> = - x p< a> (9-58) 
and 


“(01 do = z poeta (9-59) 


Equation (9-57) is also derived in Appendix I using 


the conservation equations 


It should be noted that there is no need for any 
arbitrary parameters » and 8. This is because the above 


Swmeamical equations can be determined exactly. 
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With Eqs. (9-58) and (9-59), the dynamic Eqs. (9-54) 


and (9-55) and the growth Eqs. (9-51) and (9-52) become, 
nee oe = eto 


~~ <t7> = - (9-60) 


ee ee (9-61) 





R! o 
d.<t2> + ‘iia 2c oe 2 
RR | | ont? Rie ee 





Re 
2 2 
= a2 oe (9-62) 
a 4 2 K 
Ips (9,7) oh ot a <(a,7) > = - F pdp<t*> . (9-63) 


peo -OluULIOn Of the: Moment -Kaquations which Govern. the 


Geowtco- Ol ~vuhbe Fueruscions 

ihe ecrowthn equavions are solved in Appendix kK. 
Wie interested reader might note the differences between 
Mmm@enwerowth equations for the pressure free model and the 
radiative pressure model. The most significant difference 
mapenat the curvature constant k is in general not negligi- 
ble in the region where the matter density exceeds the 
@aadiation density. Therefore, in order to avoid difficult 
integration problems, the growth equations are solved only 


for the k equal zero universe. The advantage for doing 


ie 





this is that the results can be expressed in terms of simple 
functions. (It seems that such problems presented by a 
nonvanishing k are inherent to nonuniform models. Sachs 
and Wolfe [24], for example, are also led to consider only 
Euclidean models for the same reason.) Nevertheless, the 
results for k equal to zero should have some relevance to 
k equal to plus and minus one universes. This is because 
micwaynamics olf the universe are now only slightly modified 
Mmyecurvature effects [25] and the solutions to the growth 
Bebo Ons Ought. to therefore be approximately correct for 
euea models also. 

Werte -equel. CC. Zero, snes luctuabion tree <soluvicn. 


ier Kp and Ri? given Anse pendix late 





3a 4 
KO —_ — 
R? 
and 
a 
oeee, esr 2 
Ri? = = (9-64) 


where ag ie Lon IivetramuoniwecCiatanu Ol here: (fiom a. which 
appears in the analogous equations for the radiation filled 
universe. 


Then the growth equations become 


2 2 2S eel 2 

R Inp<t > + 5 ROB<T > 52 

2 2 

= a, Recto, cree (9-65) 
a 0 
d 
aac : 
2 \ aoe? a5) F ae 

OC) big he = fas.\ = Ip<T >. (9-00) 
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The solutions to these equations are given in Appen- 


pix KK, 
(3.1) 
epee 
0 e -1 -7/2 -6 
easier Tk au FR + DoR (9-67) 
2 2 -l/2 y Zs 
<T*> = Tk - a BOR / ae 9 Dok 3 ; (9-68 ) 


4. Special Submodels 
Again, as in the radiation model, only the following 


three models are considered. 


2 
: 2 
o-7°> =.= - a : cs FoR 7? (9-70) 
d 
(3) <t2> = 3 Dok? ee 2 DoR-® (9-71) 


The same difficulty occurs with the Fo-model as 
with the Piz model Oia Une previous SECULON«. ‘That. as. pie 
is impossible to find a nonzero Bo Wasclh i Sighs thes 
fompatible with both of the Eqs. (9-70) (a, is always greater 
than zero). Hence, only the I, and D, models are considered. 
joes tO be noted that the I ,-model eS a, medel 2orm waiter 
Mieenonuniformity increases with increasing R and the Dy- 
meael is a model for which it decreases with increasing R. 


if vEqs. (9-09 ).and (9-71) "are substituted into 


Eqs. (9-60) and (9-61) the following equations are obtained. 


it 





@eew— i luctuations ) 





0 
1 ae 
Ko + 57 B(I))R = ee 
2 
R 
Z 
Chie =o a (9-72) 
R R* 
(Dp - fluctuations) 
2 
omer ee? B(D,)R° 2 2 Ouas 
2 
R 
6 2 
- 2 B(D,)R > = - & - 2 (9-73) 
R Re 
where 
IT.a 
Es ea) 
oa NG: 
and 
Dia 
ae OG 
The solutions of Eqs. (9-72) and (9-73) turn out to 
be 


(I, = fiuctuations) 








3a BU) 
SoMa = —. ae (9-75) 
R o R 
a 
RI = -¢ - 3 B(I,)R (9-76) 
(Dp - fluctuations) 
3a 20 RGD.) 
Kp = - _ 7 + (ont) 
R R 
a 2AE CDS 
R! = aa es 32 ow a (9-78) 
R 3 R 
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ante 


Mmaeries, 


is used to 


then 
type 
sileyel 


next 


Les 


K and 


the same classification scheme, which is depicted 
and 2, is used to classify these solutions that 
classify the solutions for the radiation model, 
egeidyvotound that: the Ty - fluctuations are of 


ae Do - fluctuations are of type M. The discus= 


of this model and the radiation model is given in the 


@imapver . 


ie 





X. SUMMARY AND CONCLUSIONS 


fee VALIDATION OF PROCEDURE 

Pe, Ore discussing the results of this study recarding 
mee development of a cosmological model to demonstrate the 
applicability of a moment hierarchy scheme to perturbation 
Beebiems of general relativity, it would be well to show 
that the solutions obtained are compatible with ones obtained 
byea completely different procedure. 

felimwmown ian the Uiterature 15 the work by Sachs and 
Wolfe [24] concerning the anisotropy of the microwave 3°K 
Meexeoround radiation due to perturbations of geometry. In 
tors Study a scheme is presented for analyzing small, linear 
Merturbations away from an unperturbed, k equal zero 
Robertson-Walker metric. The method of solving for the 
memcurbations is to fives take the spatial Fourier transform 
eee perturbations, solve them in kK Space, “and then. trans— 
mame toem back into position space. They do this for both 
Mee cust and radiation models. 

For the dust model (pressure free) they find two kinds 
of terms. In both types, 6p decreases although 69/p increases 
for one and decreases for the other (in their notation, o 
is the mass density and 6p is its small perturbation). 
Furthermore, the characteristic time for their relative 
mgereasing solution is of the same order as the character-— 


[merc time for the background metric. They report their 
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solutions for 69 to be of the form 


6p = pR7?/2 + ie (10-1) 


where D and I are functions only of the spatial coordinates. 


3 


Because 9 goes as I/R-, 


38 = pR7 3/4 LoTe (10-2) 


Mere the constant of proportionality is absorbed into D 
ema l. 

As shown by Eq. (9-57), the ratio of é6p/p must be propor- 
lmomal tot. Hence, according to Eq. (10-2), t must also 


have the form 


he ie (10-3) 
hams, t’ must be given by 

apa hea oe oe (10-4) 
Bopeeq. (9-68) says that 

-—OF 
yy -3 0 .-l/2 2 
2 a a So ne — 
auc > 5 Dy R + 5 R + Ik (10-5) 


He 21S seen that the exponents of each term agree with 
those of Sachs and Wolfe. This agreement is noteworthy 
because Eqs. (10-4) and (10-5) are produced by two very 
Omeeterent procedures. 

ier determination of the characteristic time T-. for the 


i 


Iy (increasing) perturbation term t can be obtained from 


Z 
mie ratio of <t*> to S00 8G) 2, ona is. 


8) 





: aL as 





7 s ae eee (10—o) 
T 2 
e- <(9,t) > 
where c is the speed of light. Using Eq. (9-69) 
2 R3 
Ty - (10-7) 
Caan 


The characteristic FAme Tp for the background geometry can 


be found from the eae OsOL R? to ORE That is, 





2 Le 
a SS (10-38) 
c2 RI? 
ere 6d. (o= 64) )i5 is becomes 
: | 
T,’ = = (10-9) 
cas 


Indeed, the ratio sp/p which 1s proportional to 7 is pre- 
dicted by the model to take place on the same time Scale 
as the packground because the ratio of Ty to Tp is one. 
This constitutes a second point of agreement Witheoe@cns 
and Wolfe. 

The above comparisons are made only for the dust models. 
But a comparison between radiation models can be made in a 
similar manner. Sachs and Wolfe indicate fhay vmey neglect 
transport processes intel derivation. Honeces ene: CCle =e 
ponding model of this “Oe tseune Bere 1eakage" model POL 
which A and B are both equal to one. 

For their radiation model, sachs ana Wolfe report that 


one pR7e (VO 20) 
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so that 


: 2-4 
aa = DER (ioe 


where D is a function only of the spatial coordinates. 
The general solution for <t*> given by Ea. (9-36) with 


moma Bp €qual to one is 


<tT*> = 


: D 
\ =f 
US tt ee Se (10-12) 


The Pa/3 solution is omitted from the general solution of 
Eq. (10-12) because, as stated in the text, certain sign 
problems are associated with it. 


But Eq. (9-24) indicates that 


<C> 


: Seats Guo eels 


Mis indicates that 


2 


5 a T2 (10-15) 
moat is, 

ay (Taya) 4 , [Puzsl 

ec |e i (10-16) 


The comparison between Eqs. (10-11) and (10-16) indicates 
that although there is disagreement concerning the existence 
of an RI Crh. (iets le ae vecieny sconce nae = ic ene Cerin. 
itis not possible to assess the reason for this discrepency 
owing to the omission of the solution procedure from Sach's 
ema Wolfe's paper. However, it is significant to this work 


mae With all the possibilities for the values of the 
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exponents due to the many possibilities for )} and g, there 
isean agreement concerning the existence of an Bee term. 

Thus, there is a major agreement between the functional 
Mmomms Of both the k equal to zero, dust and radiation filled 
models. Because this agreement is accomplished by two very 
Samterent perturbation techniques, support is given to the 
method of this work for treating perturbations by a hierar- 
ey Ol moment equations. 

It should again be mentioned that the agreement of the 
radiation model is contingent upon the selection of the 
"Zero leakage" values for A and 8 because Sachs and Wolfe 
memme that the velocities are zero. It is reasonable to 
eee 110W Valid is this assumption. In the derivation of the 
meeel, the attempt was made to include the effects of the 
memeocities by the-introduction of the parameters ) and @ 
whose magnitudes are estimated in Appendix H to be of 
Omer one. 

Figures 3 and 4 give the types of solutions that can be 
egeceted for all possible combinations of A and gp. In the 


mere OL the I = Type peruurvar ion, the model classifica- 


ie 


imo for A and @ values of order one is clearly type Kk. But 


mene Lhe D - type, both the K and M - type solutions are 


Wes 
meeciple, in fact, the zero leakage values of A and 6 are 
observed to lie on the boundary between the K and M ~ type 
eemmcions. Thus, if a Days - type model is being considered, 


@meneglect of the velocities is not trivial for very differ-— 


eme dynamical situations could occur even if the velocities 
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Bieeewextremely small. But for an 1473 = UVOe . bie eter Lee Goi 
the velocities does not change the solution type. In this 
case the approximation of zero velocities seems to be justi- 
ied. indeed, Sachs and Wolfe consider only this type. 
mewever, A and ® not only affect the classification but 


also determine the exponent (which is j+4 for D and it4 


Ls 
ior Ty 73) OF Rein ithe expression or <1->.)) 1f Pigmre 5 
tg@eswlines of constant i and j are given in A-8 space. From 


Maese Curves, it appears that for A’ and g of order one, 


Sec iets 6. 4 iO =a) 


In the case of X and 8 equal to exactly one, j+4 is equal to 
~4 and it+t4 equal to 4. ieee estimates are conservative. 
Perhaps a more reasonable estimate would be to assume that 
Gme likelihood of A and B both being on the large side of 
gear order one magnitude is small. Then the upper limit 
might be pushed down to perhaps -6 and 6 respectively. 

mans 16 Seems that the zero velocity assumption can 
reasonably be questioned on the grounds that the exponent of 
<tT*> could be different enough from the zero velocity expo- 
Meme that a many order of magnitude change of R could produce 
4a value of <t*> that is much larger (or smaller) than the 
fem esponding zero velocity value. This can be seen more 
Mmiemcivatively as follows. Consider, for example, a radia- 
tion model with an 14/3 - type perturbation. At some value 


of R equal to R, 
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4 
rE 
et?>4 = = coe (oa) 


Pierne the sub and superscript i refers to the i-value of the 


exponent. If it is assumed that the ratio Maga is 


0 
one at R, then this ratio is in general different at other 
meenes Of R. That is, this ratio is a measure of the effect 
of a given combination of A and 8 upon the nonuniformity 


foecerms of the "zero velocity" nonuniformity. This ratio 


momcasily found to be, for an arbitrary R 


aT Goh i 
— = (10-19) 





where i = ~ 4 + 2/Y2(\tB) . 

Thus, unless A and 8B are such that i is near enough to 
Meee so that <t*> is of order ae the eftecto ot tne 
feepemsi1on through many orders of magnitude upon the nonuni- 
formity compared to the zero velocity nonuniformity is signif- 
icantly increased as seen from Eq. (10-19). If, however, 

R changes at most by only a few orders of magnitude ina 
memwer problem, the discrepency is not as serious and may be 
Memerable if i is close enough to 0. Hence, although the 


Zero velocity approximation is extremely convenient mathe- 


macacally, its validity does not appear to be on solid footing. 


B. DESCRIPTION AND COMPARISON OF MODELS 

Moe staced in SectionesA ol This chapter, there. are cyo 
mm@erinsically different types of nonuniformities associated 
fae, both the dust and radiation filled models. One iS an 
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increasing (I-type) and the other is a decreasing (D-type) 
Monuniformity. If the ratio of <t*> to <(a,t) > is computed 
for all four models (two for each aq value) in the same way 
as done previously so that a characteristic time might be 
Eseamated for the relative perturbations, then it is found 
meee, TOr all models, this characteristic time is of the 
Same order as the characteristic time for the average 
background. 

Because of the existence of the two D and I-type pertur- 
Mmarons, it is natural to ask which kind characterizes our 
eaeverse. This question can only be answered through obser- 
feeton. Although such observations are difficult to perform, 
thewe seems to be some suggestion that the I-type of pertur- 
bation is the one persian to Our” Own Universe, ~~ One 
ifeason why this might be true is that the average background 
metric mero De Mexpanding more: raprdily Chan the supercilusver 
ipemenich our galaxy belongs. <A second reason is that, in the 
case of D-type perturbations, galaxies might be discarded 
from their otherwise clustered configurations [24]. 

eieancing the i-type dust model, it is possible to estimate 
the parameter To Which eoverns the amplitude of the metric 
perturbation <t*>. This can be accomplished by using the 
data reported by Partridge and Wilkinson [19] on the aniso- 
tropy of the so-called cosmic blackbody radiation. They 


3 


report that the ratio 6T/T is of the order 10 ~ where T is 
memes temperature of the radiation (3°K) and 6T is the fluctua- 


tion of T with angle of observation. Sachs and Wolfe [24] 
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give a very simple expression that relates ¢§T/T to bp ,/P, 


where p, is the current mass density. 


O0 
ST eee 2 
Ww 7 HO (gh) = 
ane 
where H, = y— (' denotes derivative wrt. x°), R. is the 


0 0 
current value of R, and L is the characteristic length of the 


Peeturbations. By Eq. (9-57), Eq. (10-20) can be written as 


2 
2 ee ee <2) (10-21) 


9 alee 
(H,L) 


The constant I. can be estimated by means of Eqs. (9-69) 


and (10-21). That is, 





7 (10-22) 


tH 
R 
M 
£ 
a 
O> 


a 46) 1 
where 
etn ie ; (10-23) 


ne HpR is taken to be one and if the factor 16/9 is neglect- 
0 


eapesance it is of order one, Eq. (10-22) becomes simply, 





OS R 2 
0 


we 6T)* 2 
I {| (10-24) 


If the fluctuations are assumed to be the size of about 
moO MPC (24, 31] (one MPC is equal to 3.26 light years) 
meeenat ¢ is of order ten (R, ~ 10% em), then for Partridge 


and Wilkinson's estimate of 6T/T, Eq. (10-24) gives for Ly 
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oe 
Ip = 10 aon | (10-25) 


Using this value for Ty; it is-interesting to-deternine 
the effect a term of this strength would have upon the average 
dynamics if it is assumed that the long wavelengths are 
muoortant to the expansion characteristics [35]. It is well 
known that the k = 0, unperturbed uniform universe model 
expands forever [2]. (This can be seen easily by integrating 
Eq. (9-64) to find Ra t°/2.) But owing to the additional 
mourece term, the dynamics of the model is perturbed. As 
Pommnved Out earlier, the solution for the k = 0, Io model 
fo tLype K. Because this means that R is perturbed below 
the uniform solution, (see Fig. 1), the question naturally 
perces whether or not there €xists a maximum value for R. 
me answer this, Ri must be set to zero in Eq. (9-76) and a 


Positive value of R greater than R, must be determined. ef 


mis 15 done, such a value ie is found such that 


1/2 
Soe a (10-26) 
0 
Using Eq. (10-25), 
R a nO CM, (10-27) 


maxX 


That is, if the observed anisotropy of the blackbody radia- 
mien is due to perturbations of space as large as super-— 
elluster distances and are of the Ig - type, then the 


"effective energy" of these longwave fluctuations is 
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sufficient to prevent an unlimited expansion of a k = Q, 
@vst tilled universe which would otherwise expand forever. 
Tt should be pointed out that this model is not valid 
mer all values of R. For small R, the radiation density, 
which goes as ae overwheims the mass density, which goes 
as R7?, peas ipelmt at 1S me. Honeer Valiageyo wepearnel a 
aust model and it iS proper to use the radiation model for 
smaller R-values [15]. On the other hand, for large R, 
<t*> approaches one. As it does so, the assumption of small 
amplitude fluctuations which was used to truncate the moment 
hierarchy, becomes invalid. The point at which this happens 
can be easily determined. With <t?> equal to one, Eq. (9-69) 
gives 


al 


_ ie 
(I) 


(10-28) 


where Ro is that value of R for which <t*> is equals vo one: 
iirc ( 10-20 )) tsecompared to Eq.a.)(l0=26) "ie is=feund 

that a WATChrdnedteavese that sta bhes, vurning pedo.) 

the fluctuations are as large or larger than the average 

metric. However, it should be noted that because the model 

iemnot applicable at the turning point itself, it may nov 

maerefore be justifiably used to determine nae ot LAS 20 at ie 


emlty is not serious though since the ratio of the “extrapo- 


" , = R : 
lated value" of R_., given by Ea oe 26) to , is 


: h 
Max = ~ (10-29) 
io Se 


30 





Pico) 1 S501 order one. Tnus the real value of Raa should 
still be of the order of 10% ecm as given by Eq. (10-27) 
since the theory breaks down close to the estimated value 

of ees Se It should finally be noted that if other models 
are considered which also have “turning points," (either as 
a minimum or maximum) the weak fluctuation assumption breaks 
meu just beiore the turning point occurs. Thus, if the 
miayvsics Of an assumed turning is to be accurately discussed, 


emourong fluctuation model should be used. 
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APPENDIX A 


The Reynold Stresses 


The Navier-Stokes equation for an incompressible fluid 


ese | 2 | 
au > = i a 
ee a = VP + yy2V (A-1) 
= 


mecre V and P are the velocity and pressure fields relative 
to a fixed coordinate system and v is the ratio n/p where 
Meas the viscosity coefficient and p is the mass density. 
meieincompressible flow, po is a constant throughout all 
ace . 

imeorager tO pain some insight into the perturbation 
Mmeenmique that is presented in this paper, the reader should 
look upon the Neealere coe: Squat On as 711 so were. timc 
Einstein Field equations. The left hand side of Eq. (A-1) 
mmo oecous FO the geometrical side of the field equations 
memerne right hand side is analogous to the mass side 
(that is, the side which contains the energy momentum ten- 
mem). The velocity field v plays the same role in the 
Navier-Stokes equation as does the metric tensor Sy aie ele 
Eemstein field equations. 

im fluid turbulence theory, the following decompositions 
are assumed [4]. (In this appendix., the equation after the 


semicolan refers to the analogous equation in the text.) 


>. 
Veze<V> + Vv (A~2); (4-9) 


oe 





Peace (A-3); (4-9) 
where < > denotes "average." 
If Eqs. (A-2) and (A-3) are substituted into Eq. (A-1), 


the following equation results 


ae ae > > 
=> 
< 1 — + <V> + V<V> + Vv « V<V> 


© 
“A 





=> 
= * Vepso= = Vp + v¥2<V> + vV?2¥. (A-4) 


ma fluid turbulence theory it is assumed [4] that if 


A is a random field, 





ee 
Ses =. yee (A-5);(B-2) 
and 
PL OM EHS CN=6 CBee) 
Fl i 
ox OX 


Because of Eqs. (A-2) and (A-3), 


-> 
<> 


0 (A-7); (4-14) 


and 
<p> = oa: (oo 


Then, Eq. (A-4) becomes 
=> 

acy > 

ong 





> > a 
+ <V> « V<Vo = - o Vers 


> 7 > > 
+ vVV¥S<V> - <v + Vv>. (A=9) (6-8) 306-9) 
It should be noted that Eq. (A-9) is symbolically iden- 
> -> 
meeal with Eq. (A-1) with V and P replaced with <V> and <f>. 
The only difference, however, is the appearance of the extra 


i 





source term from the "velocity side" which results from the 
nonlinear V : vv term, It 18) in: general, nonvanishing. 

This extra term is called the Reynold Stress. It can 
be interpreted as an additional source term (in addition 
bowcvhne pressure and viscous sources) which acts on the 
average flow field (<V>). Such an interpretation can be 
made because the left hand side of Eq. (A-9) is symbolically 
identical with the left hand side of Eq. (A-1). 

iieorder TOvobtTain tne functional form of the Reynaid 
mieesseS, a hierarchy of equations can be built up from the 
Navier-Stokes equation of Eq. (A-1) by multiplying it with 
mac1ous orders of v and-averaging. “BUG yUnesnontinear Vemm-: 
ee already seen, always produces a moment of one order 
meecaver than the order of the moments Which comes from the 
other linear terms. Thus,there is always one more unknown 
lan there are necessary equations available to solve the 
BEymevem Of equations. This problem is the well known closure 


meuoplem of turbulence theory which is discussed in the 


ext {3}. 
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APPENDIX B 


Definition of the Averaging Operation 


There are various definitions of average which are pos- 
Sible. Many theorists who work in the field of fluid turbu- 
iemee UuSe an integral expression for their definition of 
average. For example, Curle and Davies [4] use a time 
integral. This is an acceptable definition for experimental 
work but as pointed out by Lumley and Panofsky [13], such 
em integral expression is not always easy to define in a 
meaningful way. They suggest that an ensemble average is 
pmooanly the most basic. Thus, the notion of an ensemble 
Eyperape 1S taken as the definition for the averaging process 
[ee Or 6this model. 


B are the 


Suppose D(x',&) is a stochastic process where x 
hour space-time coordinates and € is the ensemble parameter 
warech designates Das a possible outcome for an experiment 
Meee also Ch. IV). Then <D> is defined as 

D(x"; Oe es see ig Gt D(x"; Ex? 
Qa ed (B-1) 
N>3o© N 


Mand is not dependent on 


where <D> is only a function of x 
—~. This is the same definition as that used by Lumley and 


Panofsky. 


Because this definition involves a summation, the follow- 


ing equation is assumed valid. 


oD 





doe ta. <0. De z 
\ per (2) 


SUP POCe ull onde c ied. On MiOW EOC eprOobeabwl ty 61 1 iieciemc 
D(x") falling below a certain value 6. This can be deter- 


mined by counting the proportion of members of the ensemble 


mor which D(x") = 0.) From this’a function F €anwbe-dertined 
as follows 

F(6;x) = P{p(xh) < 6} (B-3) 
[mere P stands for “Probability that." The probability 


fiemsity Lunction is defined by 


IF(63xh) 


f(D;x') = = 


(B-4) 


Then the average of any function n{D(x")] can be written as 


mice ie =| n(ejr(é,x Jad 4 (B-5) 


—OO 


96 





AL een Ac 


An Implication of Spatial Stationarity of the Random Field 


An implication of statistical spatial stationarity is 
rather easy to state because of the assumed form of the line 
emement. On a hypersurface of constant ae dx°® is zero so 
mmeu by Bas. (4=5) and (4-6) 

27,0 vd. 
Roy (x »X we) 


5 2 
ds,;.0 = - me  do* (C-1) 
ie 2 


2 
ry? 1 + i ni 





2 


an 


iieos pointed out in Chapter IV that R,, has the neaning of 


M 


mecal radius of curvature.” It is also a stochastic process 


mien means that at each point xh 


of space time it has 4a 
probability density function associated with it (see Appen- 
ox B). 

im order CO _ DrIne about a-simpiart 1¢eacion, Ol «the preposcda 
meael, consider the following assumption. Although the 


hyperspaces of constant x? 


are ¢leaerly nonuni?s erm. and. wide: 
BewvVarious expansions at each point, it is assumed that the 
me@obability of finding a g@iven Expansion at any point is 
masependent of position on that hyperspace. This require-— 
ment is then one of statistical spatial stationarity, a 
mencoot which is familiar in«fluid turbdulenee theory. This 
implies that the probsbility of measuring a given value of 


R,. between say Ee and R,,” + alge ought to be independent 


M ] M J 


yl 





OI the coordinates on this spatial hypersurface of constant 
time x° because this random function contains all the expan- 
mmom characteristics. Equivalently, the probability of 
Mmamadaine ¢t between IT and T + dT ought to also be independent 
Mmarmec Spatial coordinates. Hence, the probability density 


oiac > LON f(T;x") defined by Eq. (B-4) must be such that, 


; ft 
f(T3x°,x2) = £(Tsx°,x2 ) 


J a 


for every x* andx WHeCTreGatiN= soc ees mee nS rea roll yamean 
that [18] 
HQ) eae Cle”) (C-2) 


That is, f can only be a function of the coordinate x° and 
meer But Eq. (B-5) says that for any function h of the process 


moe 3c )> 


Ss fai CP). CL dt Fe (C-3) 


=O 


memee, <n(t)> can only depend upon and Motvaupom1 sie mmree 


iheetollows that 


OG 2 = 0 Cea) 


where j =r, 6, ¢. 
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APPENDIX D 


The Form of the Energy-Momentum Tensor 


The general form of the energy-momentum tensor for a per- 


meer luiad 1s given by [1,29] 


fy = DURUY a (Ga Se) (D-1) 


© 


where U! = dx¥/ds and D and P are the PRrOPeI hyorostcau Lc 
Wemotty and pressure. That is, these are the density and 
ieessure which would be measured by an observer moving 
instantaneously with the fluid. Tolman [29] defines a per- 
fect fluid to be a mechanical medium incapable of exerting 
meansverse stresses in the proper frame. From this defini- 
mien follows Eq. (D-1). 

Tolman [28] shows that the energy-momentum tensor for a 
disordered electromagnetic field, which is expressed in terms 
of the electric and magnetic field strengths, can be viewed 


meme operfect fluid for which 


oa ae D- (D-2) 


W|F 


However, if the pressure as aque to the random motion “or 
mercicles of cold dust, then the equation of state, which 


je well—known to kinetic theory, is given by [1] 


ik 
= — D- 
P 3 0 ( 3) 


oO Is 
po N 
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where V* is the root-mean-square random velocity of the dust 

peaxrvicles. But this velocity is much less than ec so the 

meessure, for all practical purposes, is essentially zero. 
The two models, therefore, can be written together in 


eme equation. 


P= oe aD (D=4) 
where a = 1/3 if universe is radiation dominated 
a = 0 if universe is matter dominated 


Now, the expression for the mixed energy momentum tensor 
a ie wiceaea. Jhis 18 easily found by means of -the metric 


Gensor er 


(uUMU - ot) (D-5) 


where 


Gee oe (D-6) 


In order to evaluate Ea. (D-5), U ana U, need to be 
Serermined. This can be done as follows. Consider a fluia 
element which moves along the trajectory x9 (x°), Along 


ech a trajectory, 


“Z 2 
ee ea(ar) || eve (D-7) 
(ax) 
Waeen the definition 
u = ei%e oo k (D-8 ) 
dx 
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equation (D-7) becomes 


ae 


oS = + (1 - u?) ; (D-9) 


ax? 
Notice that if the motion of the fluid element is along a 
null geodesic, u* must be equal to one. In order to keep 
the interval of such trajectories timelike (interval real) 
Memrnav tne local velocity of the Tluid element does not 


exceed the light yelocity, it is obvious that from Eq. (D=<9), 
Oise ce le oe (Day) 


iteeis natural to identify (uc) with the local velocity of 


M/2 


the fluid element at the point x4 because e ““do is the 


Bisjoical distance that the fluid element moves in a time 
I 





ax? ; Now define 
im di a 
a G2 ue ies, 
So Lonat 
dx? 
U as ics 


fimerei ore, 


Ce eee oe (D-13) 


But Eq. (D-8) can be written as 


2 2 2 
eae (Clean pee ot Cas a ere ceases | Spee 
Gx ax ax? 


im the quantities U,.> Up and ne are defined such that 


oy 





then 


iy 1S natural to identify u 


- live aise 


ax 

a rele dé 
0 dx. 

u =r siné eM*e ap 
Q ax 
2 2 2 2 

Uu o + Uso + ae 


J 


ioeceal velocity at ae Hence, 


by Ea. 


Ce ey 


ur = ye Mey 





r 
-~M/e2 
oars © 
Ur = Y = Up 
; eae 
as i sind “¢ 
(D-6), 
U = Y 
0 ; 
z M/e 
U = — Y 2 
U =- yreM/2,, 
2 6 
U = - yr siné eu, 
3 


OZ 


ao, Gue yuh 


(D-15) 


(D-16) 


COMpONnen Cl On Bu mc 


(D=17) 


(D-18) 





Then Eq. (D-5) becomes, 


ae =e ey 

0 02 
Se ern =a Tena u2) 

} r he r 
Ue eae oe le ee eae 

6 62 6 

= cebu al Fee =e) 

3 db c 2 p 
ie = - enn =- (D + . ae 

ce 2 . 

Too os - p2elhp2 = -—- (D + ae rellZe,, 

2 0 > 8 

Cc 
shied Seer Git 6 eT = - (D + Py r sin6 etre, 
e? > 

ec ee tte 

sl metas P Mens (D + >) PU,Up 
Cc 

} 3 

T, = r*sin?6 T = = (D + PY re Sano ue 
] ae r @¢ 

Tne eo ree ae et) eek ner (D-19) 

3 y 2 0° ® - 


Cc 


where D = DY 


2 
yu, tle y(1tu, - ao (D-20) 





APPENDIZ E 


Solution of Cosmological Equations with a 1/R" Source Term 


The two usual cosmological equations of general rela- 
meet y With Zero cosmological constant but with a ee 


source term are [1] 


2 
Son ues 














Ope am (E-1) 
ine R* 
and 
eee 2 
Pee ee 2 oe eee CHa) 
ee eA R? R R? 
where 
Poe (E-3) 
R 
P = acy C=) 
Pee C0 GE=5)) 
V Vv V | 
We 8 1G (E-6) 
a2 
t= dq (E-7) 
dx” 
The difference between Eqs. (E-1) and (E-2) is 
t 
K ok RK. 
— (P + + <= 2/--| . E-8) 
K(o + p,,) + ; (P Pp) = E ( 
C 
The derivative of Eq. (E-1) is 
1 
6kR! R') (R! 7 
iS Grey + p.) = ae = 6 [E] FY fs (CE 9) 
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mine substitution of Eq. (E-8) into Eq. (E-9) gives 


cepa) 


and one Sh! be _ 
(p P. as ae Cotp,) + re Ont: Cae 0) 


mEq. (E-8) is multiplied by R*, it becomes 


(Eee 


Tor 2) ae ee) ee One (E-11) 
Cc 


This is the conservation of mass equation where = and P 
are the mass density and pressure of the source respectively. 
iemis valid for any equation of state P(p) and BoD e 
buco With Eq. (B—5) , oe ie a6" py, and EQ. .( 5-4) 0 P taeqe-o 
Eq. (E-11) becomes 
' ' | ' ' 
MGs iene ea o ee) ep Gh o)) wee, (E12) 


The homogeneous equation associated with this equation is 


(pR?) + ap(R?). = 0, (E-13) 


mae solution of which is 


A 
p3 (a+) 

In order to solve Eq. (E-12), Eq. (E-14) may be substi- 
tuted into it assuming A is now a function of x®. Eq. (E-12) 
becomes 

30 +3 BO 
= ' a ; = 
ee o 3R R'p, (lta) (E-15) 
.. _v | 
With p_ = — , Eq. (E-15) becomes 
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(m-3(1 +o, )} 


{ = ~~ 
n (m-3 (1 +a ) J 


But Eq. (E-16) is easily integrable. It turns out that 


v (Rs 42)-m_ . (E-16) 


eee cyp36atl)-m + B (E-17) 
where 
m-3(1+a,) 
“" m-3ita) eae 


Finally, the substitution of Eq. (E-17) into (E-14) gives 


ae ne CV 
Poa GT ae 


inewsubstatution of Eq. (E=19) into Bat (B=) )veives en 


2 
expression for R' . 





M2 ees B rl V(1-C ) : Z 
R 3 | Taal in=2 a (E-20) 
Thus, 
a 3a ie “ 
KP = ~ “30gFI) ~ Gm co 
and 
b 
Ri? ea a 1 ~k (E-22) 
poatl pii-e 
where 
ee S (E-23) 
f = KCV (E-24 ) 
be = a (E-25) 


Tf vy = 0 (no source term), then the solutions become 


the usual cosmological solutions for p and R. For a radsav2on 
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iMevled Universe (a = 173) [20]. 





3a,, 
KO = - i -26) 
and 
2 oe 
R! = ye ee 46 (E-27) 


For a pressure free universe (a = 0), the zero source 


term solutions are 





38.4 
Ko = - (EB-28) 
R 
and 
a 
oe — a (E-29) 


It should be noted that the meaning of the constant a 
is different in the case of the radiation model (denoted 
by a.) than in the case of the pressure free model (denoted 


by aq). 
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APPENDIX F 
Elimination of the Spatial Derivatives 


from the Dynamic Equations for a Radiation Filled Universe 


The terms involving the spatial partial derivatives of 
T can be removed from Eqs. (9-1) and (9-2) by the substitu- 


men of Eq. (9-3). The result is 


Rt 2 
Cea = TRC se as z Coen ee 3x gq2> 


R2 


2 
Ble Sn se eae 


Ro R R? 


Cpa) 


R! 2 
pe = monte cas = Sia i Ace 7 Toe 


2 
tf t 
_ 4k Po ce UR 
RR R Ke 


(F-2) 


Bue the right hand side of Eqs. (F-1) and (F-2) are not 

more xpressions that would result from substitution of the 
peerage metric into the Einstein field tensor. In order to 
bring these equations into such a form, Eq. (F-2) must first 
be multiplied by minus five and added to Eq. (F-1). If this 
Ewmeis multiplied by three eighths, then the desired form 


results. 


3 kK 2 = 3k one ° Gh 3) 





To obtain the other equation, Eq. (F-1) must be multiplied 
by minus 5/3 and Eq. (F-2) by 3. If these two equations 
are added and divided DY e12hn04 they resu vine cova ten. 1c 
formed. 
OR pad 172 oe eee 
= "5 16 ot ea on es (diet) > 


~ = a <T.> = - eo . (F-4) 


2 
The term 9, <t*> can be replaced by using Eq. (9-4). The 


final pair of equations are 


6) 3 : Re 2 
Kp + 7 K<td> ta, crete Rie ac anaes 


t 
ae re ae eee 
R R R 
2 t 
ms ae = a K<TtTd> - i Sie oer - no et 
2 
tt t 
~-¢Xti at --h- ft ; CE aio) 
Re Re R R 
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APPFNDIX G 
Klimination of the Spatial Derivatives 


from the Dynamic Equations for a Dust Filled Universe 


Hie cl Mipimneat ion Of the terms: which compauumpne. coacial 


partial derivatives in Eqs. (9-48), (9-49) and (9-50) 


gives 
' 2 
-3xp -3R 9 <rto + p< (a,t) > - 2.«<ta@ 
Z 
+ 3k <t*> = lek = ene. + alae CG) 
R2 R2 R? sat 
Sa ae One — ae o> d g2cq?> + 2 <(9 t) > 
6 ook 0 er 9 y 
) 2 
? 1 
= 2 K<td> + pee aes a Hk + os + 3k” (G-2) 
R2 R2 Ba R 


Now, if Eq. (G-2) is multiplied by minus five and added to 
Eq. (G-1) and the resulting equation multiplied by three 


euenths, the following equation is obtained. 


GR 2 Dees 9 as 
Oe erp mas Opn e eet qe eee glen 
2 
t 
+ 2 K<Td> - p< <V- See . C= 30) 


If Eq. (G-1) is multiplied by minus five thirds and 
ba. (G-2) by three, and if these two resulting equations 


are added and divided by eight, the following equation is 


obtained. 


LAO 








tt ! 
: r a er th fe oh age EN aan 
R 


Now, the term 9*<1*>that appears in Eqs. (G-3) and (G-4) 


ean be eliminated by using Eq. (9-51). The result is 


2 ’ 
Kp + i acres > 3 RY Rui + ?, K<Td> 


{ 
z i Ne Nee US, pee ees (Ges) 
R? R? Re 
{ 
-~F<Q t) > te ait? + =e x<td> 
lk so. k OR" RI 
ieee SS ae) ee ee eee Be (G-6) 
Re Re R R 
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APPENDIX H 
The Dynamics of the Energy Momentum Tensor 


for Drsorderead Radiacien 


The dynamics of the disordered radiation fluid is 


described by Eq. (9-10) 


cae: vy = 0 


where 


av 


mHY . = UV + V HO 4 u cS 
> V a ae eee CH 1) 


dae components of ieee are Piven In Append’ x]. andthe cor 
Mouents Of pie are given by Eq. (5-10). Because of the 
summation which appears in Eq. (H-1) with respect to the 
indices y, Eq. (H-1) represents a set of four differential 
eamations. Adler [1] points out that the equation with yu 
aqua LO - Zero 1S the €nergy—-conservation equation while 

mae remaining three are the relativistic equivalent of “New-— 


mom's Law of dynamics applied to an ideal fluid. The con- 


Servavion of energy is therefore expressed by 
yates seo oh eee hg (H=2) 
Vv Va av 


With the appropriate substitutions from Appendix D and 
Eq. (5-10), Eq. (H-2) becomes 


AD (baete Shine aoe pan) ee moe enn2 | 
0 a2 2 0 02 


-M/2 E 
+e 9{((D + — )u. J (H-3) 
ec? COn ue. 


ae 





+ e209 i + =)£(D + —)u J GOnit« 
oe (H~3) 
a-M/2 1 


v 
P 94 L(D a ae 


Sa 


ae cot 8 R 
e CSF og) ir ere De oy Ug J 


~M72 1 P 
qa al 


P _ 
Coane he et ee 


Cc 


if the equation of state given by Eq. (5-22) is used, 


Eq. (H-3) becomes, 


as 29 MD 7 


Wr 


2 a 2 
9,(Dut) + 3 (a, M)Du 


y -~M/2 —~M/2 } 
+ 3 Le Gane i? 


ap Dues) + e 


~M/2 1 


COGad 
( 9.M + = )(Du,) 


a e Ve 2 a, (Du, ) + e 6 
r 


£ 


=i 
2 Tee 1 


+ een 
r sin 6 


3, (Duy) 


pliZe at 


+ — 
r sin 8 


( dM )t Dud I = Or (H-4) 


PPVMOrGer. to, la SUltadLe EXDressi Ons (lhe mi eirus replace 
Mme correlation functions <td> and <td,d> in Eqs. (9-4) and 
(9-5), the following procedure is used. 

If Eq. (H-4) is multiplied by a special function 


h(t,9,7) and averaged it becomes 


<ho, D> + 2<h(d M)D> + = <ha,(Du~ )> 


i 
: 


a de 





ee _M/2 h 
+ a <hnDu (9,M)> + U<e eae 3 (Du, )> 
hDu 
een oe 


i snoulad be noted that in Eq. (H-5) only partial deriva-— 
tives with respect to x°® and @ appear explicitly. Those 
terms which contain partial derivatives with respect tor 
midst are buried in the terms containing the ¢ partial deri- 
vatives. To see how this can be done, consider the follow- 
ing argument. In Eq. (H-4), there are three types of terms 
mavOolVvVing Spatial partial derivative with respect to r, 6, 
maa @. On a spatial hypersurface of constant x° it is 
mem olanle CO assume that the averace of a quantity associa— 
MmeGawloth a given direction in space ought to be equal To 
the average of that same quantity associated with a different 
direction. Since the r, 6, and ¢ terms of Eq. (H-4) have 
miemuregl roles. (i.e, transport ef mass in’r;190, 6 direct rons 
as Adier [1] points out) then the average of these terms 
(multiplied by the function h which has no explicit spatial 
dependence) should be equal. 

Now, D can be decomposed into its average and fluctuat- 


imeancomponents. That is, 
Dev elGras eas CRO emtclng iS) (H-6) 


pich that <D> = 9 . 
if} Ko. (H-5) 1s e@xpréessed in terms of t instead of M) 


it becomes, using Eq. (H-6) 
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ec A ee SUBIR rh NI ene : OU, 


it should 
do not create 
ifeon fluid to 


Monat is, 
Weed 
whenever 


Too 


fwhis makes certain moments of Eq. 


ecners. 


Consider now two cases. 


ewo special functions, t and 0,1. 


G San. 7 20 r sin 6 o 6 





hu = 
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< a acta < a ee > 
Ee Poi ie ae oso 2 U4 
nh rat 
Ne ae eer et 
eS ie eed et 
hdu Coal 
-M/2 g ® 2 a 
ise" psin 6 > ; eee 
be expected that weak fluctuations of geometry 


pressure gradients which accelerate the radia- 


local velocities near the speed of J1ipghc: 


CH=) 


(H-7) small compared to 


Lee ‘hisymbolicaliystand -10r 


Then 


ED: 





~ 


eS (H-9) 


Thus, if the higher order moments of Eq. (H-7) are neglected 


by virtue of the smallness of u’*,this equation becomes 





a 1 _ a 
Sol Sct “a <hd> + 2p<hd,t> + 5 o<hd,u7> 
or 2 — -M/e2 h 
+ = 0< 0.T> + ——_—__—__—__ 
3 one a Ho<e SSGERA 16 eae 
hu hu 
-M/e2 d = —, -M/2 db 
< ee te sate Se toes 
+ Axe Yr sin: 6 ee + 4p<e r sin 86 Fa 
=O. (H-10) 


As pointed out by Harrison and others [9], any perturba- 
fron of the radiation fluid characterized by d, U,,» Ugs Uy 
feel propagate with the sound velocity of the fluid. For 


fmawration this velocity is near the speed of light. 


V = (Las 


) 
sound  Yv dD Vea . 


fenee, such an accoustical signal would propagate along 
Biter ceedesic lines: That is, ds? (accoustical Signal) 


meenear enough to zero that 


dx® x A aaae = eve Ps ito dh Se. (H-11) 


Ete chis implies that the three moments which contain the 


spatial partial derivatives with respect to @ in Eq. (H-10) 


Sueht to be of the same order as those for which Bale x 
— poms Treplaccd. Dy sd.u. FOr texanip le 
isin 6 0 > 
-72 h 
a i ecingor ag og 





Then, Eq. (H-10) becomes after discarding the insignificant 


moments 





<ha d> eee: 20<ha ,t> 
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We h 


+ Unce = = EE 
Won<e Sane 3444? Or. (H-12) 


[i Srme Lirst two terms of Eq. (H=12)arerconsolidaued 


into one term f<had>(for an appropriate f,), it becomes 


M/e h 


fi<ho d> a 4on<e r sino ene 


= - 2o<ha.t> : (He43) 


it the characteristic time for the fluctuation is comparable 


mo vtvhe characteristic time - a Of the everace eeoneury 
c R 2 


Mgemimvne SeCCona Germ 15 2equalk “in magnivuce toyrne Tirse. 


Mnmtinis Case, f., which is a function of x°., is of order 


1 ? 


- 


Oe . 
Novercons der Une tremeimime equagtons TY. ys Oforey 
eueial GO One. CWO. and three A Ailer-some rauner tedious 


ealculations, these equations become 


(Te = 0) 
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state is used 
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Ferms are discarded and if the equation of 


these equations become simply 


et 


WR _ 1 -M/2 
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(H-17) 
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ee ae + <a DN + I+ ‘Du, ) 
= _ 1 eM*2 ana a oe : (H-17) 


In order to learn something about the behavior of the 


-M/e2 h 
sine 9 Uy? which appears in Ea. (H-13), 


moment <e 
Pensider the average of the third of Eq. (H-17) multiplied 
ay hh. 

ib a7 M/e h 


p<ho » Uy? = Tor anu (H-18) 


(To obtain Eq. (H-18), the higher order moments are 
merrecced.) , 

Because: Of the previous arguments ,concerningwuneare la 
Mmmonshiip Of Spatial and temporal partial derivatives, stne 


mollowing equation should hold for some f Of OTrdersone. 


— -M/2 _h _ = ; 
o<e r sine oan te ; Pi<hg,d> g (H-19) 


If Eq. (H-19) is substituted into Eq. (H-13), the follow- 


ime, equation results. 


ned = ear, poets , (H-20) 


If the changes of ad and t occur over the same time, 


me must be true that 


ade = f p<ht> Geo) 


where a is a weak function of x® of order one (unless by 


Sience tat ot eeomalis a 
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But in Eqs. (9-4) and (9-5), the correlations <td> 
and <(o T)d> are required. Since the ratio of D to D is 
(ie 73° ) (see Eq. (9-6)) whieh 1s very close to one, 
it can be assumed that. LOr  SOmesruner Lonms dp and Bn» both 


weak functions of x® and of order one as described above, 


foe 1Ollowing equations hold. 


TC 2) .p<t*> 
<(d,T)d>= - BP d<T>- (He2e) 
iy 


r and Bp are expanded in a Taylor series and only 
the constant terms are retained because of the supposed 
weak dependence of these functions on x°, the following 


epproximate equations result. 

<td> = - 2Ap<t*> 

“(0 ees Bont e (H-23) 
These constants are then to be viewed as arbitrary parameters 


of the model. Eqs. (H-23) are therefore consistent with the 


memservVation equations. 
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APPENDIX I 
The Dynamics of the Energy Momentum Tensor 


for a Matter Filled Universe 


The dynamics of the energy-momentum distribution is 


eaven by 


ee = Oe (i=) 


mine conservation of energy equation corresponding to u 
equal to zero is given by Eq. (H-3) with a, and P equal to 


zero. (It should be noted that u, and P equal to zero with 


J 
Meunspecified are solutions to the three momentum equations 
Mor which p = 1, 2, 3. See Eqs. (H-14), (H-15), and (H-16). 


maat is, 
> D + 3(,M)D = ja (1-2) 


372M 


Memeiis €quation is multiplied by e and integrated, 


meme Lollowing equation results. 





D a BoN72 ; (I-3) 
With 

Deacon O (I-4) 
and 

eo/eM Z o3/eG (141)37¢ . (1-5) 
equation (I-3) becomes for small tT 

eee = zs 
ae de (I-6) 


ee 





hele e, 





ase A 
e 
and 
poe OT (I-8) 
2 ¢ 


nat is, for a pressure free universe, the stochastic process 


mesa LUNCECLOn of therstochastic prosesset.. Hence; 


Udo. = p<t > (I-9) 


Nofw 


and 


eo le - : pont ; Cin.) 


123 





er i La) 
Solution of the Growth Equations 


Lor. a Radtarpion we led wuna verse 


2 
27>. 2 = R4 <(0oTt) > 


a5 
if yo= <T >, 5 Vee Se 1 ge UNC Er Olen 


r 
Eqs. (9-34) and (9-35) are 
2 
ee ary. 4 X se = O)VU==ex 2 (J-1) 
ax ie 

and 

dz y _ 48 ay 

ax i Xx a x" ax saz) 


+ 
ferme NR 1S a constant with the dimension of length. All 


mee variables (y, z, x) in Eqs. (J-1) and (J=-2) are 
dimensionless. 

The solution of these two coupled differential equations 
meeceeds as follows. If Eq. (J-2) is multiplied by x", 
the left hand side becomes a derivative of the function x’‘z. 


Mas equation may then be integrated directly to give 
y = 1 x'z + 7 | ie) 
4g g 


where ie oro arti LLrary CONS lanl. Ol “iit Cota emo. 

From Eq. (J-3), the first and second derivatives of y 
mewn respect to x may be found in terms of 4 and x. If 
these are substituted into Eq. (J-1), a second order differ- 


Eimgetal equation in 2 results. / 
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2 S2RAL 
ce one + 9x a a oe een = ae (J-4) 
a+ “ x" 


Equation (J-4) can be transformed into a second order 
differential equation with constant coefficients by the 


change of variable 





ae ee (J-5) 
eu S , 

a*z az _ —Ut 

ae? + 8 ae tO 2 ee ee a) gear : ES 


The associated homogeneous equation of Eq. (j-6) is 


©, 





eg dz 2 
> ee Pe ee 


ites complementary solution Za is 


- 


z, = Aev’ + Bed” Ges) 
where ese a OOK DB) 
=e alee ee 


and A and 8 are arbitrary constants. 

The complementary solution Za Cam ie used Toe finde ine 
Beameral solution to.Eq.. (J-6). This m@eudone by assuming 
the general solution is given by Eq. (J-8) but with A and 
Perunepiogns of tC. 

if “nis SOlut lon ut Orm. 18) substi tated into ba. (d—o 


pmoaitferential equation in A and B results. 
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d+ Come. 
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= s2pfe7 (J-9) 


iieceiirst tyo terms are’ zero Dy Virtue wor {tne tacu. thet 


oa and edt are solutions to Eq. (J-7). In the equation. 


Ge Ai tye." Bite a (J-10) 


immece Cependent variables of t are present: Zy A. and Bb. 
mwetn general, three variables can be made to satisiy three 
Pemavions. However, so far, only two such equations have 
been specified, noe, Eq. (J-4) and Eq. (J-10). Hence, 


Pmeiraeequation can be freely chosen. A convenwent choice 


tS 
Pe clea cade) © i 
Ae + @€ at 0 Cie ki 
so that Eq. (J-9) becomes 
att Ae 
Gey Gh de on ee cs Hit . 
ee et a 328Af ¢ (J-12) 


These two equations can be written in matrix form. 


~~ ow ~ 


He . (J-13) 
- ert ed 
where aL 
Aen ae 
ae at 
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~ ae 
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aB 
dt 
7 0 
pl 32p,r ie" | | . , 


By Cramer's Rule, 
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Now, Eqs. (J-15) and (J-16) may be integrated to give 


ee -(444)t 2 = 
” Gi aety © ue — 
32,ef oe 
Br- ong e = 


where D and E are arbitrary constants. Hence, by Eq. (J-10) 


and by the identity x = a 
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Ze ee ee x (J-19) 
Bf 
where KF = cee and PS a new arbitrary constanv. by ska. 
(J-3) 3 
Veda te D a + ee i (J-20) 


ii etnese soluvions Por Yoana Z arewsubst itucvedrinvo tec. 
(J-1) and (J-2), they are found to satisfy these equations. 
if y and z of Eqs. (J-19) and (J-20) are translated back 


2 
into <t*> and <(9,T) >, Eqs. (J-19) and (J-20) become 


Se “1/3 Rit 4 13 RI*4 (j2209 

So =~ WAP) /gR + ae + Da gh" Clee) 
Mere i = — 4 + 2V72(0\+B) 

j= - 4 ~ 2vEC0FB) 
and Ty 3: Dy /3: and Pa 73 are arbitrary constants. (14 73 = 
Dyptith Aaa Da /3 Z Eyptity 
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APPENDIX K 


solution of the Growth Equations for a Matter Filled Universe 


3 
ee ee are gone Bes: ri 
ee ah Se on <(9,1) >; x = ey (for some R*), the 
growth Eqs. (9-65) and (9-66) are 
ee te Ea: = cy ee 2 (K-1) 
: Be 
De Es 
7, eee 2 aed : ie) 


If Eq. (K-1) is differentiated and Eq. (K-2) is substi- 
Peeeead into it, a third order differential equation for 


results. 
Dee Ue Dy OO Me Gee ee eee cere (K-3) 


If a change of variables is made such that 


x =e Ce) 
then Eq. (K-3) becomes 


ara + 21Z,, + 61z, + 4Oz =0. (K-5) 


The general solution of this equation is simply, 


ae Toe a ee Kop 


In terms of x, Eq. (K-6) becomes 


1 Ay 172 


Z = Coan 5 oo 


+ Che? 


mi this equation is substituted into Eq. (K-2), an equation 


form y' “results. if this equation is integrated, the 
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hollowing expression for y- appears. 


y = Gx - oo eof oe : icc ee (K-8) 


But Eqs. (K-7) and (K-8) are general solutions to the 
set of the two coupled Eqs. (K-3) and (K-2). Four integra- 
Moms are required to solve this system which yields four 
arbitrary constants (G, H, J, and K above). However, the 
given system of equations (Eq. (K-11) and (K-2)) require only 
three integrations. Because Eq. (K-3) is derived from these 
equations by differentiation, Eqs. (K-7) and (K-8) must be 
meee eenecral solution to the piven equations With one of (tne 
eopropriate arbitrary constants equal to Zero. " By anspecvtion: 


ioe equation 


y= K and z = 0 


methough a solution to Eqs. (K-3) and (K-2), is not a solu- 
mirern tO the original equations. /Hences the @eneralesoluvicn 


to Eqs. (K-1) and (K-2) must be 


pe ee ao +: pee (K-9) 
y = Gx - = x le + = Tx" (K-10) 
Hence, ; 
Sl oet) > 7 y : 
: -~~Rt+pRl/2 4p R? CR 
a 0 0 0 
d 
and 
oF, 
<t2> = I,R° - —~ RIE + = D,R | (K-12) 
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APPENDIX L 


Slassit ication of Solutions for the Radiation Filled Universe 


2 
Withechne<exceptionm of therinverrationm Ofma. 6) Sac... -co, 
and (E-21) of Appendix E give the solutions to Eqs. (E-1) 
and (E-2). These solutions are, with k neglected as usual 


(see Eqs. (9-41) to (9-45)), 


Kp = - Shag eng ll (L-1) 
R* A 
a. -2 mate 
sey ak + b,R (L-2 ) 
where | , 
CCL 
rs = te LOX S46 = a (E=3) 
iy 3 Y2(A+B) 
he 3 TU tte ae ee ace | (L-4) 
3 ¥2(AX+8) 
3 fot: CCT) 73? 2a + V2 48) (L-5) 
14/3 36 YO(X+B) 
b 2 Even oe eo 
ry 36 2(A+8) 
m =f =-— 4 + 2Y2(0)t8B) (L-7) 
I 
73 
m Se oc OU Gus) 
Ve 


ine £ subscript designates the type of fluctuation. (15 72 - 


type corresponds to a situation where the nonuniformity 


ot 





increases with increasing R and Dy 3 - type corresponds to 
a situation where the nonuniformity decreases with inceasing 
R.) 

adhere are four kinds of solutions based upon the manner 
in which kp and R are displaced from the fluctuation free 
molutions «p* and R*. That is, f, and Dy can in principle 
Be Cach positive or negative depending on the choice of 4 
and B. 

In order to determine how a given (A,8) displace xp and 
Reirom Kp* and R*, it is necessary to locate those values 
(1,8) for which f and b are equal to zero. Then a set (\,8) 
meet tes Ona particular side of the=tUwoucurves of “2ero 
rh and “Zero Da determines the sien.or fh and Dy and hence 
the direction of displacement of kp and R from xp* and R%*. 

inorder to find the line of zero 1, IneX-Besvpace. cen 


A 
mecer first the following definition. 


Wie oy 2 OGG Jikan (L~9) 
The set (\,8) for which f, equals zero is given by the 
equation 
28 
Oye ee ee (L-10) 


~ f2(X4+B) 


maere the plus sign is used if A represents Dy and the 
meus sien if A represents 11/3" This equation as it stands 
memecocthner Gifftficult to plot directly. »Bul Une plotting is 
simplified if Eq. (L-10) is separated into two parametric 


equations with parameter n. 


se 





The two parametric equations 8@(n) and iA(n) are found to 


be 


(L-11) 


WD 
fl 
hops 


5n* - 
Sees 


bo} po 
eee” 


vg ae) 2 ai 
8 fe 


= 
n | gn? Fn — 1 | (1-12) 


[i +h 


where n 2 0. The top sign corresponds to the D3 —- type 
fluctuation and the bottom sign to the 1473 - type 
mime tuat ion. 

By varying n, the equation 8(A) can be traced out. This 


eeievion i1seplotted in Pigs. 3 and 4 for both if and D 


1/3 1/3 


type fluctuations respectively. 


URCVOCUAtLTON-OL US VSI88 71 Or ype = UTS MuCh sip eer mre 


A 
mea.  Wailth straightforward alecebra;, thas equation is 


B = A4(2d - 1) | (L-13) 


Hor both the 14/3 and Ds /3 —Jtype fluctuations —~ <this ava 
tron is also plotted in Figs. 3 and 4. it should be noted 
meat Bq. (L-13) is quadratic in X. Hence, according to this 
Polaron. “pnere Should be two ValUucs sor ©) -correspondine co 
eeeciven Value of @. However, only one value Satisties the 
mequirement that by equal Zero. I[havetbwere varcotrvorya iesren 
X comes from the fact that Eq. (L-13) is derived by a squar- 


mm operation which thereby introduces an extra root. 


io 





APPENDIX M 


index and Location of Symbols Used in Dissertation 


thes o lowing list.comcains hems ymbols- ccd“ Chis 


work. 


ReEcompanying each symbol 1s the location-where tne 


Metaniclon or explanation of the symbol may be found. 
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Figure 1. Classification Designations 
for the Increasing Type Fluctuations. 
mages sOlid lines are the fluctuation free solutions for R 
and op. The broken lines depict the possible deviations of 
meena «€6©fplh6ufrom these unperturbed solutions. The letters K, 
meer, and N label the four possible combinations of perturbed 


Pelucions for R and o from their unperturbed values. 
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PIgure 2. iClassifieation Decicnacions 
for-the Decreasing Type: Ff lucuuations:. 
Mae solid lines are the fluctuation free solutions for R 
mmecieO. lhe broken lines depict the possible deviations of 
need 0 from these unperturbed solutions. The letters K,; lL, 
fend N label the four possible combinations of perturbed 


fPoemutions for R and op from their unperturbed values. 
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Rapure 3. Solution Types in A-8 Space for the Increasing 


Ev peur Mictwaclt Ons for the Radiavaen en iiied  Uniwer sc 


mne 4-8 plane is sectioned into four regions. The perturbed 


solutions for each region for R and oe are depicted as broken 


The fluctuation free solutions are represented as 


solid lines. 
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Figure 4, Solution Types in A-8 Space for the Decreasing 


Myoe “iucvuataons for whe Radiation filled Universe. 


mige A-B plane is sectioned into three regions. The perturbed 


solutions for each region for R and po are depicted as broken 


Pnestluwctuation free solutions ares represented as 


solid lines. 
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